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ABSTRACT 

Two solutions are developed for the two dimensional 
problem of bonded linearly elastic half-planes. The first 
is for two bonded isotropic linearly elastic half-planes of 
different elastic properties having a crack along the in- 
terface as well as a perpendicular crack in one of the half- 
planes which may intersect the interface crack. The appro- 
priate integral equations are developed through the use of 
displacement dislocations in conjunction with Mellin trans- 
forms. The resulting pair of singular integral equations 
is solved by obtaining the relationships between these solu- 
tions and the weight functions for Chebyshev and Jacobi 
polynomials. The second solution is for two bonded iso- 
tropic linearly elastic half-planes of different elastic 
properties having a crack along the interface, as well as a 
perpendicular crack in each of the half-planes, either or 
both of which may intersect the interface crack. The ap- 
propriate integral equations are again developed through 
the use of displacement dislocations in conjunction with 
Mellin transforms. The resulting three singular integral 
equations are solved by obtaining the relationships between 
these solutions and the weight functions for Chebyshev and 
Jacobi polynomials in a similar manner to the pair of 
equations in the first solution. 

For each solution, numerical results are presented for 


Ill 


the stress intensity factors ^ strain energy release rate, 
stresses and displacements. 

The behavior predicted by the above studies was in- 
vestigated experimentally using polymers for the material 
pairs. Very close agreement was found for the critical 
stress intensity factor at fracture for the perpendicular 
crack near the interface. Fracture along the interface 
proved to be inconclusive due to difficulties in obtaining a 
brittle bond. Some interesting and predictable behavior 
regarding the potential for the crack to cross the interface 
was observed and is discussed. 
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Introduction 

In considering the fracture of composite materials due 
to the presence of imperfections, one finds a common mode of 
failure to be that of a crack originating in one material 
then extending until it reaches a bonded interface and spread 
ing along the bond line. Other possible directions of crack 
growth are for the crack to cross the bond line or to be 
reflected by the interface into the original material. 

An investigation of the effects of a crack in the near 
vicinity of a material interface is then essential in the 
fracture analysis of composite materials. In the analytical 
study of this problem a specific geometry that leads to 
tractable integral equations is the case of two elastic half- 
planes containing a finite length crack in one of the half- 
planes with the crack being perpendicular to the material 
interface. Of particular interest is the behavior of the 
stresses as one end of the crack approaches the interface 
and either terminates, crosses into the adjacent material, 
extends along the bond without crossing the interface, or 
upon reaching the interface, extends into the adjacent ma- 
terial and debonds along the interface. The instances of 
the crack terminating at the interface Or crossing the inter- 
face have been investigated by Erdogan, Cook, and Biricikoglu 
in [1] and [2]. The present study considers the possibility 
of the crack spreading along the interface. 
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The appropriate integral equations are developed through 
the use of displacement dislocations in conjunction with 
Mellin transforms. The resulting singular integral equations 
are solved by obtaining the relationships between these solu- 
tions and the weight functions for Chebyshev and Jacobi 
polynomials. 


Formulation 

Assuming a coordinate system having the origin at the 
center of the interface crack, as shown in Figure 1, it is 
seen to be convenient to represent the solution in terms of 
the polar variables (r, q ) , as the crack surfaces, the in- 
terface, and the horizontal plane of symmetry lie on constant 
6 surfaces. 

An effective means of formulating linear elasticity 
problems having boundary conditions of this type is to make 
use of Mellin transforms on the radial variable. It is 
demonstrated in [1] and [2] that the use of displacement 
dislocations in conjunction with Mellin transforms gives a 
very simple and straightforward technique to develop the 
integral equations for a large class of problems of this 
general form. 

The half-planes are assumed to be loaded with uniform 
stresses a , a and a as shown in Figure 1, with the 
stresses being related in such a manner as to give constant 
strains in both the x and y directions at points remote 
from the cracks. The following relations must then hold. 
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0=0 
1 2 


and 


0=0 
1 2 



for generalized plane stress 


(I - 

(1- V^)E^ 


VgCH- 
(1 - V^) E 

\ X '' 2 



for plain 
strain. 


The complete solution may be represented as 

a total = Oj + Oj-j 

where = stresses in the half planes without cracks and 

loaded at infinity 

and o = stresses in the half planes having no applied 
II 

loads at infinity and applied stresses on the 

crack surfaces equal to the negative of those 

given by o . For the loads as stated above 
I 

this would require normal stresses of and 

on the interface and perpendicular crack 
respectively. 

The solution for a is developed below in terms of 

II 

a general system of applied tractions on the crack surfaces, 
although still requiring symmetry about the y = 0 plane, 
and will be restricted to the above constant normal stresses 

I 

only for the numerical results. 

As the present study is closely related to [1] and [2], 
an attempt to follow the form and notation of those investi- 
gations will be made, where possible, in-order to eliminate 
unnecessary duplication. 
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Following [1], the Mellin transforms of the stresses 
and displacements are: 


M[r^aj^(r,e)] = 2j^(sj0) = 2 i(s+l)X 


is0 


A^se +Bj^(s+l)e 


i(s+2)0 _ -i(s+2)0 


-®k ® 


M[r®Tj^j.j.(r,e)] = - s(s+ 1) 


is0 _ -is0 
A^e + A^e 


- (s+l)(s+^) 


i(s+2)0 _ -i(s+2)0 


Bke + B^e 


M[r^Vk,(r,0)/2[Xk ] = Vk(s,0) 


s+1 


is0 

A^se + Bk(s+l)e +Kt, Bi,e 


( 1 ) 

i(s+2)0 ~i(s+2)0 

B o 

k "k 


^k = 


where 

t 

3— 4v^ for plane strain 

(3“Vk)/ (l+Vj^) for generalized plane stress 
*> 

and |j,k j are the shear modulus and Poisson’s ratio 
respectively. 

The transformed functions above are [1] 

' = ■'krst’''®) + ^■'kee(''’®> ’ 

■'krr(r,e) = 

''k(r>e) = 2l"kL ^ ^ = 


and 


( 2 ) 
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for region 1 and 2 respectively. 


Using displacement dislocations'*' on the surfaces 0 = ^ 
and Q = 7X , the boundary conditions are: 


T o(rjO) = u o(rjo) ,= 0 j on 0 = 0 in material 2, (3) 

2^0 2c7 


= 0 

^ u^g(r,7T) = 


, on 0 = 7T in material 1 


- l/2h(r)6(r-r ) 


(^) 


and 




^ u^^(r,7T/2)- ^ u2j.(r,7T/2)= -f(r)6(r-r^) 

^ u^g(r,7T/2)- ^ u^^(r,7T/2)= -g(r)6 (r-r^)^ 

By defining the unknown functions , k = Ij 2, 

A=C+iC , B=C+iC,A=C+iC.B=C 

1 1 21 3 4 2 5 6 2 7 

and on transforming the above boundary conditions 
suiting eight equations specifying the unknowns 
as follows: 


interface 

e = 7 t/2 . . 

( 5 ) 

as 

+ iC 

8 

j the.re- 
are 


^ Displacement dislocation as used in the present text 
implies a step discontinuity in the displacement slope at 
a particular point. 
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sC 

+ (s 

+ l)c + 

0 

11 

0 


(6) 



6 


8 

8 





sC 

+ (s 

+ 1)C - 

• K C = 0 ^ 


(7) 



6 

8 

2 8 


sC^ sin(TTS 

) + 

sC 

2 

cos 

(tts) + ( 

s + 2)C^ sin (tts) 








+ (s + 2)C cos(tts) = 0 
4 

(8) 

sC sin(TTS 

) + 

sC 

cos 

(tts) + ( 

s+ 1 - K )c sin 

(tts) 


1 

2 


1 3 

s+1 





+ 

(s + 1 - 

K )C COs(tTs)= ■ 

9 

(9) 





1 4 

2(s+l) 


sC cos (tts/2) -sC sin( 7 Ts/ 2 ) -sC cos (tts/2 ) + sC sin ( 7 rs/ 2 ) 

1 2 3 4 

- s(C -C )cos(7Ts/2) = 0 f (lO) 

5 7 

sC sin (7Ts/2) + sC cos (tts/2) -(s+ 2)C sin( 7 Ts/ 2 ) 

X 2 3 

-(s+ 2)C^ cos(Trs/2)-[sC^ - (s+ 2)C^] sin(7Ts/2) = 0 , (ll) 


- m[ sC^ cos (tts/2 )- sC^ sin (tts/2 ) - (s + 1 + )C^ cos(tts/2) 

+ (s+l+K^)C^ sin (tts/2)] 

u,pf (r«) 

+ [sC^-(s+ 1+ k^)C_^]cos(tts/2)= - — / (12) 


- m[sC sin (tts/2 ) + sC cos (tts/2 )- (s + 1 - k )C sin (tts/2) 

X 2 X 3 

-• (s+ 1 - K )c COs(7Ts/2)1 
M-2g(r )r®‘*'^ 

+ [sC -(s + 1 - K )C ]sin(TTs/2)= , ( 13 ) 

^ 27 s + 1 


where iti = M-^/m-^ and in Equations (lO) through ( 13 ) use has 

been made of Equations ( 6 ) and ( 7 ) which require that C =C =0. 

6 8 


The remaining functions are found to be 
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g-|- 

cos(tts) (s+ 2)h(r^)r^ 

C = - C + — 

- ^ - - - ’ - -sin (tts ) , 2s (s -^- l) (1+ K^)sin(7Ts) 


(1^) 


s+1 

LL„r 
^2 o 


( s + 1 ) ( 1 + mK ) (m + K ) 


sin (tts/ 2)[(s + l)(m+K^) 


~ ( 1 + ItlK ) ] 


+ g(r^)cos (tts/2) [ (s+ l) (m+ k:^)+ ( 1+ mK 
S+ 1 


+ 




2s (s + l) (l + K ) 




(1-m) (s+ 1) (2s+ 3) (mK^-Kg) 


(1+ mK ) 


+ 


(m+ K ) 
a 


] 


(15) 


^ , s + 1 

cos (tts) 

C = - C 

^ sin(TTs) 2 (s+1)(1+k )sin(7Ts) 


( 16 ) 


s+1 

C = 

^ (s+l)(l+mK ) 


+ 


[f (r^)sin(7rs/2) + g(r^)cos(7Ts/2)] 

^^i(l-ro) (2s +3)h(r^)r^ 

2(s + l)(l+ K^)(l+ mK^) 


(17) 


C = 


U 3^ 

^2 O 


S+1 


s(s+l)(l+mK ) (m + K ) sin (tts) 


[■f (ro)sin(7Ts/2) X 


[(s+ 1) (1+ mK^)-(m+ K^)]-g (r^) cos (7 Ts/2) X 


[ (s + 1) (l + mK^)+(m + K^)] J + 


^2Mr,)r;-^" 
2s (s + l) sin (tts) 


i2s+,3). _ (s + 1) 


(1+mKj^) (m+K) 


( 18 ) 
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and 


s+ 1 


C = 


O, 


(s + l) (m + K ) sin(Trs) 


[ f ( r^) sin (tts/2 ) -g (r^) cos (tts/2 ) ] 


•^2 'o' O 


2(s + l)(m+K )sin(7Ts) 


(19) 


Substituting into the first of Equations (l), the transformed 
stresses on the surfaces 0 = tt/ 2 and 0 = tt are 

s -j- 1 

11„(1+ mK + m +3C )r f 

Z (sjTt/ 2) = ^ f(r) tan (tts/2) 

^ (1+ mKj(m+K^) L 

- ig(rQ)cot(7Ts/2)+7 [g(r^) -if (r^) ] 






2(1+ kJ 


r (l-m)( 2 s+ 3 ) (mK -K ) 1 

i 2 (s+ 1 ) 

L (1+ mK^) (m + K ) J 


(m + K ) 
2 


(l-^)(2s+3) (mK^-K^) 


[ 2 ( 3 + 2) + 

sin(7Ts/2) L (l+mK^) (m +Kg) Jcos(tts/2) 


1 — 

J cos (tT! 


and 


( 20 ) 


s +1 

ILL r 

S (s,tt) = 

( 1 + mK^) (m + Kg) 


f(^o) X 


[2s(m+ K2)-(1+ mK^)+(m+ K )] + g(r ) X 

^ cos (tts/2) ° 


[ 2 s (m+ K^)+ ( 1 + mK^)+ (m+ K )] 


sin (tts/2) 


+ 


2iHih(r^)r^ 


s +1 


(1 +K^) 


cot (tts) - 

- (mK^-Kg) 

(1-m) (2s+ 1) (2s+ 3)3 

1 

+ 

- 2 (m+ K ) 

\ 2 ' 

2 ( 1 + mK^ ) 

sin (tts ) 


( 21 ) 
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where in Equation (20) 

(m+ K )-(!+ mx ) 

y ^ 

(m + K )+ (1+ itiK^ ) 

and in Equation (21) 2 ^(s,7t) = M[ ir^T^^gg (r, r^,7r) ] j as 

The integral equations in terms of the unknown density- 
functions f(r), g(r)j and h(r) are given by the conditions 
that 


and 


T(r)+iN(r) = r 

o 

( 22 ) 


N^(r) = rT^gg(r,r^,7T)dr^ ( 23 ) 

o 

where N(r) and T(r) are the normal and shear stresses on 
the interface and N^(r) is the normal stress on the sym- 
metry line. 

On noting that 


and 


f(rQ) = g(r^) = 0 , c < r^ < ~ 
h(r ) = 0,0<r <a,b<r^< 


O' ” ^ ' ~o ' ^ ~ ' “O 

then inverting Equations (20) and (21 ) and substituting 

into Equations (22) and (23) one has 
u,^[ 1+ mx, + m+ K^l r , 

N(r)-lT(r) = — — | - i (5)+ig(5)] ^ 

(1+ mx^) (m+ x^) L ^ o ^ ^ 

+ /^[f (O-ig(C)] - 7 [f (t) + ig(r)]J 
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+ — 

(1+ K^jTT a 


( 


2(t|^ - r^)- 


(1-m) - 3r^) j ^ 


1 + IllK 


. ItlK, -K„ , 

^ + — 1 + ir 2(n^ - r^)- 


2 2 
ri + r m + K 


1 "'^1-^2 


(1-m) (3 ti^ - r^) 

( 1 + ) 




dn 


2 , 2 m + K_ 

T) + r 2 


2 2 
•n + r 


(24) 


and 


N (s) = — ^ — / Ti(ti) X 

a ^ 2(1+ mK^) 


r 8s^ 

12s 1 

^ p 3(l-m)“ 

1 

L (ti + s) 

(n + s)^J 

”2 m + K ’*’1+ mK, 

1-2 

T) + S 


2m-. 


+ 


/ 

o 


.tt( 1+ mK^ ) (m+ K ) 


-sf(|) (l + mK^) + 


(iti+Kg) (5|^+ s^)-’ 

.2 2 

I + S 


r ^+K 2 ) (31^- s^)T 

+ |g(?) 1^(1+ n»K^)- — 2 J 


d| 




. (25) 


Equations (24) and (25) are valid for all values of r and 
in particular 

N(r) - iT(r),= p(r)-iq(r) , 0 < r < c , (26) 

and 

^o^®^ ^ Po^®^ ^ a < s < b , ( 27 ) 

where p(r)^ q(r) are the applied normal and shear stresses 
on the interface crack and Pq(s) is the applied normal 
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stress on the perpendicular crack. Therefore, with the 
above restrictions on r and s in Equations (24) and (25), 
one has the integral equations necessary for the specifica- 
tion of the unknown density functions in terms of the known 
applied stresses. Note, if h(t) = 0 , Equation (24) with 
0 <( r <( c is the appropriate integral equation for an in- 
terface crack [3] and similarly if f(t) = g(t) = 0 with 
a <( s <( b , Equation (25) is the integral equation for a 
crack perpendicular to a bonded interface, [1, Equation 
(4.7)]. Equations (24) and (25) are similar to [2, Equa- 
tions (7.a,b)]. However, in the present study no difficulty 
arises in letting a = 0 as was the case in [2], That is, 
the second integrals in Equations (24) and ( 25 ) do not con- 
tribute to the singular behavior of the density functions 
as long as c ^ 0 . This is seen by noting that if a = 0 , 
following [4], then h(t) = H*(t)/(t-b)^ , where H*(t) 
satisfies a Holder condition on the closed interval 0 <( t <(b. 
The nature of the singular stress field near the crack tips 
is then the same as in the individual problems and the 
effect of the two cracks is only to change the value of the 
stress intensity factors. The proof of this follows from 
[4, p. 75], as 

f(t)+ig(t) = 0(t) = (28) 

(t-c)*^ 


and 
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h(t) = 


H(t) 
H^(t) 


for a ^ 0 


(29) 


for a = 0 


(t-b)^ 

with G(t), H(t) and H*(t) satisfying a Holder condition 
on the appropriate limits. Then letting 


$ 


(^) 1 ^ dt = f dt 

' 2-1X2. ■' t-z 2 -ixx -L /. „\ 

o o (t-c;^(t-zj 


"h 


G(c) 


X r ^ 

^ o (t-c)^(t-z) (r-c) ^sin (ttI ) 

+ Bounded function at r = c , (30) 


dt 


Q(g) 


i7T^ 


and 


near z = c , with z = r , one has 

’ G(c)cot(7T|) 

<£)(r) = - ^ ^ — + Bounded terns 


2i(r-c)'" 

Substituting into Equation (24)^ and multiplying by (r-c)^ 
with r-»-c , one finds the equation for | to be 


(31) 




Similarly, for a )> 0 , let 

1 .b 


Hill 


h^z) ^ f 

^ ^ ■ 2vri i (t-a)“(t-b)^t-z) ^ 


(32) 


(33) 


from which 


h(z) = 


H(a)e 


iTia 


H(b) 


2i (z-a)°^(b-a)^sin(7ra) 2i (b-a)*^(z-b) ^sin (ttP) 

+ bounded terms . (34) 

Substituting into Equation (24) and multiplying successively 
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by (z-a)*^ , z - 

are 

and 

or^ therefore^ 


then (z-b) , z b the two equations 

cot(Tra) = 0 , 

cot(TrP) = 0 
a = p = 1/2 . 


If a= 0 , h(t) = and p=l/ 2 j a= 0 in Equation ( 34 ). 

(t-b)P 

Making the following changes in Equations ( 24 ) and (25): 

a/0 


^ = ct, r=cy, a^ = (b-a)/2/b^= (bfa)/2, t) = 


o o 
bt , a= 0 j 


a x + , a / 0 

s = ^ P(y)=p(cy), Q(y) = q(cy), 

bx , a = 0 , 

<5i> (t) = g(ct) + if (ct) ^ ^(t)=h(Tj) , P^(x)=p^(s)^ and noting that 


[even function in t] dt = - / [g(t)+if(t)] 


o 

and 


-c 


[even function in t]dt 
1 c 

f^tg(t) [even function in t] dt = ■^ / t[g(t)+if (t) ] )( 

o -c 

[even function in t]dt, 

the pair of singular integral equations are normalized as 

1 


0 (t) - 70(y)= - ^ [Q(y)-iP(y)]+ I 0(t) c(t,y)dt, 

- y < 1 ( 35 ) 


and 


1 ^ d-t- ^ ^ 

- f f + f 0 (t)H(t,x)dt = -/ 0 (t)l(t,x)dt 

-1 -1 T . .. 


1+ K 


+ o., ^ - 1 < X < 1. (36) 


2|lj^ O 

The functions G(t,y) ^ H(t,x) and l(tjx) are 
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G(t,y)= - 




+ cy 2(ri^-c^y^)- 


a3_7T(l + Ki) 

2(n^-c=y^; 

(1-m) (STi^-c^y^ri 1 
(l+iriK^) 


r * 

Ff/p P p^ (l-m)(n^-3cV)-| 1 mKi-Kg"! 

/iT) 2(n -c y2) ^ ^ __+ ^ 

[_L (l+mK^^) j n +c^y^ nH-Kg J 




nH-K 


2 J 


ri^+c^y^ 


(37) 


a* 

1-m 

r 8s^ 

12s 1 

pmK^-Kg 3(l-m)‘j 

1 

2tt 

1+mKj^ 

L(rrfs)^ 


L nvfKg 1 + mit^ J 

Tj+s” 


(38) 


I(t,x)= 


cm(l+K^) 


is |^(l+mK^) 


(m+K ) (5c^t^+s^) 


+ 


27T (l+mic^ ) (itvfKg) 

(nH-K.) (3c^t^-s^ 


5c t^+s*^) 'I 

t^+ J 


+ ct 


(l+mK ) 

L 1 c^t 


where a = 


(l+ItlK^ + ItH-Kg) 
1 “ ^2 (1+mK^) (nvfKg) 


2+ 


and a* = 


c^t^+ 


(39) 


a© if a^O 
b if a = 0 


It is of interest to note that the above equations may 
be written in terms of only two independent combinations of 
elastic constants as shown by Dundurs in [6], The following 
definitions prove convenient in the present work. Let 

k = 

then^ 


1 + mK^ ni + 1^2 

^ - 1 ■ V ^ = T~rz~ 

1 1+K 2 1 + K_ 


Hi 


kik^ 


a 


(1 + K^) (k^+ k^) (k^ + k^ -1) 



l6 


l-m 

1+ mK^ “ ^ 

niK^ - Kg - 1 

m + Kg " kg ^ 


m(l+ K^) kj^+ kg -1 

m+ K„ “ k„ 

2 2 

m(l+ K^) + ^2 

1+ itiK^ ~ k^ 


and 

ki- kg 

^ = mr 

^ X 2 

The usefulness in making such a change is most evident 
if one is able to write an asymptotic expansion for the solu- 
tion in the vicinity of the crack tip. For example see 
Ashbaugh [7] for an analysis of the problem considered in 
[ 1] . A similar investigation is being undertaken for the 
present problem. 

Equation ( 35 ) is a singular integral equation of the 
second kind, and is similar to [8, Equation (l4)]j although 
somewhat simpler due to the manner in which the unknown 
function, 0(t), appears. Following the procedure of [8], 

[9], [10], and [11], it is possible to determine 0(t) in 
terms of the unknown function . f{t) and the loads. Sub- 
stituting into Equation ( 36 ), for 0(t), one then has a 
singular integral equation of the first kind of the same 
form as [11, Equation (4, .11)]. For a = 0 Equation ( 36 ) is 
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an equation of the second kind as the singularity at the 
origin vanishes. The solution will first^be developed for 
a )« 0 and then a = 0 . It should be noted that the de- 
termination of 0(t) is the same in either case; that is, 
the first part of the following analysis is true for all 
values of a . 

Referring to [8] one can write 


and 


0(y) = w(y)^(y), w(y) = (l-y)“(l+y)^ , I Y I < 1 

1 1 1 1+"Y 

a=-2-io5, p=-- 2 +ico with oi = -^ . 


Noting, as in [8] , that w(y) is the weight function of the 

(a, a) 

Jacobi polynomials ^ solution is then written as 


00 (a^a) 

0(y) = w(y) Z 


n=0 


( 40 ) 


with the Jacobi polynomials satisfying the orthogonality 
relation [8] , 

1 ■ (ot.a) (a, a) 

/ w(y)P^(y) Pjjj(y) dy 

-1 

= 0 , n ^ m , 


(a+a+1) 

2 .r(n+g+l)r 
(2n+ a+ a + l)nl r 


= cosMtioJ ’ n = m = 0 . (41) 


|n + a^ pVl) = n = ">?^0 , 


The displacements are then found by integrating Equa- 
tions ( 4 ) and (5), from which the constant and the con- 

tinuity condition of the function ^(t) are specified. 
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Integration of Equation (4) gives the opening displacement 
-on 0 = 7T as- follows: where if a = 0/ a^ is replaced by b/ 

= - -f- f{t)dt + -f- I ^(t)dt . (42) 

~i -1 

Similarly, integrating Equation (5) and separating real and 


imaginary parts gives, 

y 

u (yj7r/2)- u (y,7r/2)= - dm / 0(t)dt + Im [C*] , 

^ ^ o 

and (43) 

(y.7r/2)- u^g(y,7r/2)= - cRe 0(t)dt + Re [C*] , 

o 

1 

where C* = c/ </)(t)dt . (44) 

o 

It follows from the form of 0(t) that 

Im [C^] = - ^ /%(t)dt , (45) 

-1 

then as 

u^j.(o,7r/2) - u^^(o,7t/2) = u^^(o,7t/ 2) = - U^g(-1,7T) , (46) 


or, therefore, from Equations (42), (43), and (46), 

u (o,7t/2) = - ^ I 0(t)dt = - -|- / 0(t)dt . (4?) 

-1 -1 

Now using Equation (40) and (4l) the constant is 

ib cosh(7ro)) i 

^o = ■ 7TC /^0(t)dt , (48) 

where Equation (42) requires that 

1 

/ 0(t)dt =0 if a. ^ 0 . ( 49 ) 

-1 

The opening displacement, u (o,Tr/2) = - u (-1,7 t), is 

ir 10 

then given in terms of as 
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u,^(o,^/2) = - 2 coilM 


(50) 


where is seen to be imaginary from Equations (42) and 
(48). 

Now substitute Equation (40) into Equation (35) j =>-nd 
using the relation [ 8 , Equation (22)] and [10] with n ^ 0 , 


(a,p) 


(a.P) (-QC.-P) 


^ /^w(t)P^(t) - 7 w(y)P^(y) = -^§4- Pj^_^(y) , I y I < 1 

(51) 

one hasj for n )> 1 , 


(-a,-p) 


^fi" ^ ^ [Q(y)-P(y)] + /%(t)G(t,y)dt. (52) 

n=i 


-1 


-a -6 

Multiplying Equation (52) by (l-y) (l+y) * then 

integrating the resulting expression with respect to y 
over [- 1 , 1 ] for k = 0 , 1 , 2 ,..., and using the orthogonali- 
ty relations. Equation (4l), the complex constants are 

then given explicitly in terms of integrals of the unknown 
function ^(t), with given by Equation (48), as 


C = 


2i 


^l- 7 ^L(n-l,-a,-p) 


I 

-1 


- [Q(y)-iP(y)] +/^V'(t)G(t,y)d^X 


(l-y) “ (l+y)"^ n=l,2,... 

(53) 

if ■^(t) = 0 , which corresponds to an interface crack only, 
and letting P(y)= -a^,Q(y) = 0, Equations (49) and ( 53 ) give 

2 a^ 




(-a,-p) 


Co= 0 , Cn= 0 , n >2 and 


C, = 




(54) 


1- 1' 
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This is the same as [8^ Equation ( 26 )] except for the al- 
gebraic sign which is due to the difference in sigh of 0 (t) “ 
i.e.j see Equation (5) above and [8^ Equation ( 4 )]. One 
similarly finds the stress intensity factor, as in [8], to 
be 

k^+ik^= lim ^/c (y-l) ““ (y+l) . 

y--l+ 

From Equation ( 35 ) for y -*• l"*" and noting that the integral 
on V^(t) is bounded and </)(y) =0, lyl )> 1 then 


Q(y)-iP(y) = T - iT = - cti / 0(t) (55) 

ire lee ^ _i t-y 


or as 


°° ™ n (a,p) 00 -iTra r) 


then 


(56) 


1 1 ^+- -iTTa ^ Q 00 

;^/<^(t):^=e (l+ 7 )(y-l)“(y+l)P 2 C P (y) . 

-1 ^ n=0 


(a.P) 


n=0 


2a, 


For 


^1 = 


a, 


' =o = °' '^n = ° " > 2 


iflfl+i^ire " °^i ^-7^ (y-i)“(y+i)^ c^pjy) 


(a.p) 


■99 


and as 


(a,P) 

Pi(y) = 9 (y- 2 ico) , 


then + ik^= V^a^(l-2ioo) , (57) 

which is the same as [8, Equation (31)]. Note misprint in 

[8] for k^ . In general the stress intensity factors are 

^ 00 (a.P) 

k^ + ik^= 2 

n=0 


(58) 
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Now substituting for 0 (t)j using Equations ( 40 ) ^ ( 48 ) j 
and (53) j into Equation (36) one has^ 

^ + / ^(t)["H(t,x)- cosh(7T03) / w(z)l(z,x)dz 


-1 


Pi 00 

+ — 2 


-1 L. 
1 


Vl_^2 n=l L(n-l,-a,-p) -1 


1 (a.P) 

/ w(z)l(z,x)P^(z) dz X 


/ G(t,y)(l-y) ( 1 +y) ^ dyj dt = 


+ 


2i 


ay 1-7 


00 -1 JL (*^ J P) 

- 2 / w(z)l(z,x)P (z) dzX 

2 L(n-l,-a,-p) -1 


/jQ(y) -iP(y)](l-y)"“(l+y)'^ 


(-« 3 -p) 

If Q(y) = Oj P(y)= - ^ and as Pq(7) =1 with 


( 59 ) 


k(x,t) = H(tjx)- cosh(TrcD)/ w(z)l(z,x)dz 


Pi 00 

+ ^ 2 


CTT 

1 


-1 


Vl-72 n=l L(n-l,-a,-p) -1 


1 (ct,P) 

/ w(z)l(z,x)P^(z) dz X 


I G(t,y)(l-y) “(1+y) ^ ^n-i 


-P 


(-a,-p) 


(60) 


then Equation ( 59 ) is 


-1 


n ^ ^ V'(t)k(x,t)dt = Pq(x) 

2a X ^ (a,p) 

/ w(z)l(z,x)p^ (z) dz 


-1 


a 


y 1-7^ -X 


(61) 


-ib ^ 

^ The term cosh( 7 To))/ w(z)l(z,x)dz is only present if 

a =0 , if a ^ 0 this term is set equal to zero. 
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which is the same form as [11, Equation (4.1)] where the 
right-hand side is denoted by g(x) . 


For a )> 0 this is an equation of the first kind and the 
solution follows exactly as in [H], in terms of the 

Chebyshev polynomials. Let 

1 1 

^(t) = F(t) with / Tj/{t)dt = 0 , ( 62 ) 

-1 

and further assume Pq(^) = “ • 

Substituting into Equation ( 6 l) and referring to [11], the 
solution is 


N 

S 

k=l 


N [t;^ + W ' ’^k)] = - “i 

20_ 1 (a, 6 ) 

- ; / v/(z)l( 2 ,x )P (z) dz, 

-1 ^ ^ 


( 63 ) 



N 

and Z g F(t, ) = 0 • . 

N is chosen large enough to give sufficient accuracy in the 
Gauss -Chebyshev integration formula, [12, p. 889 ]. The 
solution of Equation ( 63 ) then gives the value of the un- 
known function F(t) at N points on - l<(t<(l, from 
which one can now determine approximate values for from 

Equations (48) and (53) by numerically integrating the in- 
tegrals containing ^(t) . The stress intensity factors 
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for the ends a and b are then 


and 


k(a) 

k(b) 


2^3^ ^/^ 


[ 1 ] , for a )> 0 
F(-l) 

F(+l) . 


(64) 

(65) 


For a = 0 j the appropriate change of variable for rj and s 
in Equations (24) and (25) is shown in [13] to be 

■q = tb , s = xb and ^ ( t ) = h ( tb ) . 

It is further suggested in [13] that the range of definition 

of the function ^(t) be extended into the interval (-1,0) 

in order to use the Jacobian integration formula associated 

with p(“^/^-’ (t) . 
n 

Let, following [13] j 

^(t) = F(t) = — ^ F(t) ( 66 ) 


from which the solution is given by the set of equations 


N 

2 

k=l 


21SH-1 ‘ ^ '"k 


F(tj^) 
1+ K. 




r 

2 a 




- 


a Vl--Y^ 


1 X 

/ w(z)l(z,x^)P^(z) dz ( 67 ) 


where 


1,2, . . . ,N;,X^ = cos(g^^) ; ti^= cos 


ff2k- 1 
4n+2 


- 1)tt 
2 


The stress intensity factor at the end, b, is 

>'('=) = - TTir • 


(68) 


The strain energy release rate for the interface crack 
is given as [5j Equation 71]? [1^1? 

Su _ ^(1+ "■»=,) + 

Sc 2ii2[ ( 1 + mKj^)+(m+ Kg) ] ^ 2' * 


(69) 
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Niimerical Solution and Results 

The numerical solution to Equations (62) and (63) for 

a )> 0 , and Equations ( 67 ) for a = 0, has been obtained 

for the particular case of aluminum-epoxy half-planes for 

ease in comparing with the results presented in [1] and [2], 

The results presented are for constant pressure on the crack 

surfaces and no loads at infinity, or therefore the solution 

previously denoted by . The pressures are chosen by 

assuming solution to correspond to a unit vertical 

stress at infinity in the aluminum and a uniform horizontal 

stress giving zero horizontal strain in the aluminum. The 

stresses a , a and a are then, for plane strain, 
o 1 2 

Vi 

3T^v~ “ 0 *^ 2857 ^ cr^= 1.0, a^= 0.26844, for aluminum- epoxy 

(70) 

and 

V 2 

Co= = 0 . 42857 i 0*26844, 0^= 1.0, for epoxy-alumin\im 

where E szl.Oxlo”^ psi, E =4.45x10^ psi, v, = 0.30, v =0.35 

• 

With F(t) known, the constants are computed from 

Equation (49) and (53) ^nd the solution is complete. The 
displacements are then given by Equation (42), (43), (44), 
(50), the stress intensity factors by Equations (58), (64), 
(65), and the strain energy release rate by Equation (69). 

The niimber of points, N, taken in Equations (63) and 
(67) and the niimber of constants , M , computed from 

Equation (53) were taken sufficiently large to give less 
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than one percent error in any of the above computed func- 
tions. The maximum error was determined by taking succes- 
sively larger values for M and N and examining the change 
in the various computed functions mentioned above. Suffi- 
cient accuracy was assumed when the most slowly converging 
of these functions ^ which in all cases was the strain energy 
release rate, numerically converged to within the prescribed 
degree of accuracy. For example, if a= 0, b= 1, c= 1, 
values of M = N = l4 were found to be sufficient. If a = 0, 
b=l, c=0.1it was necessary to take M = 30 and N = 20. 
Convergence was somewhat more rapid for the aluminiim-epoxy 
pair than for epoxy-aluminum, apparently due to the in- 
creased distortion in the epoxy- aluminum. 

Of particular importance in the investigation is the 
potential for further extension or arrest of a crack origi- 
nating in one half-plane, as the crack either crosses the 
interface into the adjacent, material or extends along the 
interface. By assuming each possibility for extension, 
computing the appropriate stress intensity factors as a 
function of continued growth, and comparing either stress 
intensity factors for the through crack or strain energy 
release rate for the interface crack with the critical 
values, one can determine if continued extension is possible 
or if the crack will arrest after a specified growth. It 
should be noted that for non-symmetrical loading or geometry, 
possible directions for crack extension are not as restricted. 
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As an illustration of this criterion, the case of a 
crack crossing the interface first will be presented using 
the results of [1] and [2]. The stress intensity factors 
for a crack originating in material one and growing toward 
and crossing the interface into material two are presented 
in Figures 2 and 3. Figure 2 is for an aluminum-epoxy pair 
and Figure 3 for an epoxy-aliiminum pair. From Figure 2 it 
is seen that after crossing the interface into region 2 the 
stress intensity factor, K(a) is initially a decreasing 
function of continued extension and then changes slope and 
approaches the full plane solution. It then appears pos- 
sible to have some extension into the epoxy and subsequent 
arrest. For the epoxy-aluminum pair, Figure 3 indicates 
that the crack will extend to the left rather than toward 
the interface. However, if the crack does cross the inter- 
face, the stress intensity factor, K(a) grows unbounded 
with continued extension and arrest is not possible. 

Figure 4 represents the analogous cases as above, using 
the results developed in the present study, with the crack 
now assumed to extend along the interface rather than cross 
into the adjacent material. Until the crack tip, a, reaches 
the interface and begins to extend along the bond, the 
stress intensity factors, K(a) and K(b) depend only on 
and have the same fom as the left side of Figure 2 and 
Figure 3 - Therefore, in Figure 4 , only the strain energy 
release rate for the crack extending along the bond is 
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presented. Figure 4 is seen to represent the same general 
behaviour as K(a) in Figure 2 and again indicates the pos- 
sibility of arrest. It is of considerable interest to note 

^ u 

that for all material pairs investigated ^ c — ► 0 

for — = m 1 and oo as c — 0 for m 1. This 

behaviour can be anticipated from the results of [1^ Table 7] 
where the opening stress on the interface near the crack 
tip is shown to be compressive for m <( 1 and tensile for 
m > 1. 

The values of K(b) corresponding to the geometry and 
loads of Figure 4 are shown in Figure 5. The decrease in 
K(b) with increasing half length d as seen on the right hand 
side of Figure 5 is due to the closing affect of the in- 
terface crack on the 0 = tt surface. In fact^ for the 
epoxy-aluminum pair with the loads given the crack on 0 = tt 
was found to partially close at c = 2.75. The closing stress 
is shown in Figures 6 and 7 ^or aluminum- epoxy and epoxy- 
aluminum respectively. 

Some characteristic displacements for the crack surface 
on 0 = TT are depicted in Figure 8 and Figure 9 and numerical 
values for some of the significant computed functions given 
in Table 1 and Table 2. The effects of increased half-length, 
c on K(b) and the opening displacements are clearly seen in 
these two figures and the tables, with the epoxy-aluminum 
having the more noticeable change. 
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Figure 1. Bonded Elastic Half-Planes with 
Perpendicular Cracks. 
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Figure 2 


- d — ^ + 

Stress Intensity Factors For a Crack Originating 
in an Aluminum Half-Plane and Growing into an 
Adjacent Epoxy Half-Plane (from [1] and [2]). 
Plane Strain. 
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Figure 3 
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Stress Intensity Factors for a Crack Originating 
an Epoxy Half-^Plane and Growing into an Adjacent 
Aluminum Half-Plane (from [1] and [2]). Plain 
Strain. 
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Figure 5* Stress Intensity Factor, k(b), for an Intersecting 
Crack with Extension along the Interface. Plane 
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a 

c 

kM 

— 

1 

0 

1 

k 

“ - 1 

k 

“ 2 - 

STRAIN 
ENERGY 
RELEASE 
RATE 
X 10 "^ 

STRAIN 
ENERGY 
RELEASE 
RATE 
X lO"^ 
(a=b=0) 

0.3333 

0.0 

0. 6050 

0.6242 

- 

- 

- 

- 

0.111 

0.0 

0.7522 

0.8926 

- 

- 

- 

- 

0.0697 

0.0 

0.7904 

1.0210 

- 

- . 

- 

- 

0.0476 

0.0 

0.8154 

1.1377 

- 

- 

- 

- 

0.0 

0.0 

0.9582^ 

CXC*- 

- 

- 

. 

- 

0.0 

0.005 

0.97^6 

0.0 

0.0430 

0.1916 

2.3902 

0.0580 

0.0 

0.010 ■ 

0.9814 

0.0 

0.0483 

0.1979 

2.5726 

0.1159 

0.0 

0.015 

0.9861 

0 

• 

0 

0.0541 

0.1820 

2.2361 

0.1739 

0.0 

0.050 

1.0025 

0.0 

0.0873 

0.1181 

1.3379 

0.5797 

0.0 

0.075 

1.0090 

0.0 

0.1060 

0.0971 

1.2808 

0.8696 

0.0 

0.100 

1.0136 

0.0 

0.1224 

0.0829 

1.3558 

1.159^ 

0.0 

0.150 

1.0197 

0.0 

0.1510 

0.0640 

1. 6680 

1.7391 

0.0 

4.0 

0.7073 

0.0 

0.8575 

- 0.1111 

46.3612 

46.3770 


*From [ 1 ], 


Table 1 , Numerical Value for Stress Intensity Factor and 
Strain Energy Release Rate for Intersecting 
Cracks. Aluminum-Epoxy. Plane Strain^ b = l. 
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0.3333 


0.1111 


0.0476 


0.0 


0.0 


0.0 







*From [ 1 ]. 
Table 2 . 


k(b) 

a 



0.0 

0.5551 

0.5397 

0.0 

0.6109 

0.5224 

1 

0.0 

0.6172 

0.4895 

0.0 

0. 6200 

0.4605 

0.0 

0 . 6241 * 

0 . 0 * 

0.005 

0.6245 

0.0 

0.010 

0.6251 

0.0 

0.015 

0.6258 

0.0 

0.050 

0.6328 

0.0 

0.075 

0.6387 

0.0 

0.100 

0.6448 

0.0 

0.150 

0.6577 

0.0 

2.7 

0.3101 

0.0 


0.7023 



STRAIN 

STRAIN 


ENERGY 

ENERGY 

k 

RELEASE 

RELEASE 

- .2 

RATE 

RATE 


X lO”^ 

X lO"^ 
(a=b= 0 ) 






















0.1192 0.3747 9.5856 0.0580 


0.1262 0.2801 5.8517 0.1159 


0.1338 0.2035 



0.1739 


0.1571 0.0561 1.7252 0.5797 


0.1656 0.0559 1.8959 0.8696 


0.1725 0.0658 2.1136 1.1594 


0.1842 0.0804 2.5062 1.7391 





31.2066 31.3040 


Nxunerieal Values for Stress Intensity Factors and 
Strain Energy Release Rate for Intersecting Cracks. 
Epoxy -Aluminum. Plane Strain, b= 1 . 

















































































BIBLIOGRAPHY 


T. S. Cook, and F. Erdogan, "Stresses in Bonded Ma- 
terials with a Crack Perpendicular to the Inter- 
face," Int. J. Engr. Science, Vol. 10, (Nov. 1972) 

pp. 677 - 697 . 

F. Erdogan, and V. Biricikoglu, "Two Bonded Half-Planes 
with a Crack Going Through the Interface," Tech- 
nical Report, NASA-TR-72-4 , June, 1974. 

F. Erdogan, "Simultaneous Dual Integral Equations with 
Trigonometric and Bessel Kernels, " Z. Angew, Math. 
Mech., Vol. 48 ,( 1968 ), pp. 217 - 225 . 

N. I. Muskhelishvili , Singular Integral Equations . 

P. Noordhoff, Groningen, The Netherlands, (1953). 

F. Erdogan, and G. Gupta, "The Stress Analysis of 

Multi-layered Composites With a Flaw," Int. J. 
Solids Structures, Vol. 7^ (1971) ^ PP. 39-61. 

J. Dundurs, "Effect of Elastic Constants on Stress in 
a Composite Under Plane Deformation," J. of 
Composite Materials, Vol. 1, No. 3j (July, I 967 ), 
pp. 310 - 322 . 

N. E. Ashbaugh, "On the Opening of a Finite Crack 

Normal to an Interface," J. Applied Mechanics, 

Vol. 40, No. 2, (June, 1973), PP. 533-540. 

F. Erdogan, and G. D. Gupta, "Layered Composites with 
an Interface Flaw," Int. J. Solids Structures, 

Vol. 7, ( 1971 ), PP- 1089 - 1107 . 

L. N. Karpenko, "Approximate Solution of a Singular 

Integral Equation by Means of Jacobi Polynomials," 
PMM, Vol. 30 , No. 3, (i 960 ), pp. 564 - 569 . 

F. Erdogan, G. D. Gupta, and T. S. Cook, "Numerical 

Solution of Singular Integral Equations," Methods 
of Analysis and Solutions of Crack Problems . 

Edited by G. C. Sih, Noordhoff International 
Publishing, Leyden. 

F. Erdogan, and G. D. Gupta, "On the Numerical Solution 
of Singular Integral Equations," Quarterly of 
Applied Mathematics, (January, 1972), pp. 525“53^. 



40 


12. M. AbramowitZj and I. A. Stegun, "Handbook of 

^fethema^ical Functions National Bureau of 
Standards, Appl. Math. Series 554 (1964). 

13 . G. D. Gupta and F. Erdogan, "The Problem of Edge 

Cracks in an Infinite Strip," J. Applied 
Mechanics, Vol. 4l, No. 4, (December, 1974), 

pp. 1001 - 1006 . 

14. B. M. Malyshev, and R. L. Salganik, "The Strength of 

Adhesive Joints using the Theory of Fracture," 
Int. J. Fracture Mechanics, Vol. 1, (1965 )j 
pp. 114-128. 



BONDED ELASTIC HALF -PLANES WITH AN INTERFACE CRACK 
A^ A PERPEi®icuLAR Intersect iNG crack' T^ 

EXTENDS INTO THE ADJACENT MATERIAL^ 


By 

James G. Goree 

Associate Professor of Mechanics and 

Mechanical Engineering « 


Willicim A. Venezia^ 
Graduate Assistant 
Ph.D. Candidate in 
Engineering Mechanics 


The problem considered in this paper is a direct con- 
tinuation of the preceeding study to include the case of 
the crack crossing the interface as well as debonding 
along the interface. The notation of the preceeding paper ^ 
refered to as Part I, will be used throughout and citations 
from Part I will be made directly by equation or reference 
number which are followed sequentially in the present 
study. 


^ Submitted to the l4th International Congress of Theo- 
retical and Applied Mechanics. 

^ Currently at Johns Hopkins University^ Applied Physics 
Laboratory j Laurel , Maryland 20810. 



42 


^ . * .Formulation 

As Sliming a coordinate system having the origin at the 
center of the interface crack, as shown in Figure 10, .it 
is seen to be convenient to represent the solution in terms 
of the polar variables (r,0), as the crack surfaces, the 
interface, and the horizontal plane of symmetry lie on con- 
stant 0 surfaces. The cases of intersecting cracks (a or 
a* equal to zero, separately or simultaneously) and non- 
intersecting cracks (a and a* not equal to zero) are con- 
sidered. 

Following Part I, the Mellin transforms of the stresses 
and displacement are; 


" 2j^(s.e) = 2i(s+l) 


X 


r is0 i(s+2)0 _ -i(s+2)0 I 

I Aj^se + Bj^(s+l)e — J j 


-S0 -is0 


r is 

M [r^Tj^^^(r,0)J = - s(s+l) Aj^e + A^e 


- (s+l) (s+4) 

M [r^Vk(r,e)/2|ik ] = V^(s,0) 




i(s+2)0 _ -i(s+2)e 


+ 


]. 


= -(s+l)/[x. 


is0 i(s+2)0 _ -i(s+2)0 


A^^se + B^(s+l)e + Kj^Bj^e 
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where = 


for plane strain 

(3-Vj^)/(l+Vj^) for generalized plane stress 


and j are the shear modulus and Poisson's ratio 

respectively. The transformed functions above are [1]^ 


Tj^rr^^^®) = Tv^^(r,e) 


krr ' 


and 


v^(r,S) = 2H|^L + i with k = 1,; 


for region 1 and 2 respectively. 

Using displacement dislocations on the surfaces 9=0, 
9 = TT , 6 = 1^/2 , the boundary conditions are: 




^ ^i,fl(r,0)= - ~ h*(r)6(r-r^) 


■ 2 


and 


^ U^g(r,7T)= - I h(r)6(r-r^) 

T^gg(r,7T/2) = T2gg(r,7r/2) 

— U^j.(r,7r/2)- ^ U2j.(r,7r/2)= -f(r)6(r-r^) 
^ U^g(r,7r/2)- ^ U2g(r,7r/2)= -g(r)6(r-r^) 


► , on 0 = 0 in material 2 

(71) 

j on 0 = 7T in material 1 


(72) 


V interface 
f , 0 = 7 t /2 


(73) 
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where the second of Equations ( 71 ) is the only change from 
Part I. On transforming the above boundary conditions, the 
resulting eight equations specifying the unknowns cl^ are 
as follows: 


sC» + (s+l)Ct + C» = 0 , 
6 8 8 


[jLgh^r ) r 

sC« + (s+l)C» -K C» = ? ° ® 


S+1 


2 8 


( 74 ) 

2(s+ 1) 

sC» sin( 7 rs)+sC»cos( 7 rs)+(s+ 2 )c» sin( 7 rs)+(s+ 2 )C«cos( 7 Ts )=0 , ( 76 ) 

1 ■ 2 3 ' 4 

sC' sin( 7 Ts)+sC> cos (7rs)+ (s+1-k )g» sin ( 7 ts)+ (s+l-K )c cos(tts) 

IL 2 13 14 

, , V S+l 

- £(s+ 1) ’ 

sC’cos (■rrs/2)-sC* sin(-n-s/2)-sC»cos (7Ts/2)+sC« sin (tts/2) 


(77) 


- s(C* -C' ) cos (tts/2 )+s (o' -C' ) sin (tts/2 )= 0 , (78) 

'57 68 

sc' sin(7rs/2)+sC' cos (tts/2 )-(s+2)C' sin (tts/2) -(s+2)C* cos (tts/2) 

1 2 3 4 

- [sC -(s+2)c' ]sin(TTs/2)-[sC' -(s+2)C’ ]cos(tts/2)=0, (79) 

5 7 6 8 

-m[ sc’ cos (tts/2) -sC’ sin (tts/2) - (s+l+K )c’ cos (tts/2) 

1 2 13 

+ (s+l+K^)C^sin (tts/2) ]+[ sC^- (s+l+Kg)C^] cos (tts/2) 


s +1 


P-2f(r^)r„ 

-[sC*-(s+l+K )C' ] sin(TTs/2) = - “Ti , 

6 2 8 X 

-mf sc’ sin (tts/2)+sC' cos (tts/2) - (s+1-k )c' sin (tts/2) 

1 2 X 3 

- (s+l-K^)C^COs(TTs/2)]+[sC^-(s+l-K^)C^]sin(TTs/2) 


( 80 ) 




s +1 


+ [sC -(s+l-K )c ]COs(tTs/2)= -rr 

6 2 8' S+1 


O' O 


( 81 ) 


The functions c! are given by c! = C. + C* where the C. 

1 111 1 
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are listed in (l4) to (I 9 ) and C* is the change in due to 

the crack on 0 = 0 . The C* are given below .as . . . 

i 


C* = 
1 


Cl = 


c* = 

3 


c* = 

4 


c* = 

5 


c* = 

6 


c* = 

7 


c* = 

8 


C* cos 

2 Sin(TTS) 

s + 1 

2 s +3 ^ 

2 (s+l)s 1 +mxi nvt-K^ J 

C* cos ( us) 

4 sin(7Ts) 

|i^h«(r^)r^ 

2 (s+1) (l+mx^ ) 

|jL_h*(r )r 

2 (s+l)sin( 7 rs)s (l+K^) ^ 


(s+l) ( 2 s+ 3 ) ( 1 -m) 

1 + iriK, “ ra+ 

1 2 

s+l 

H 3 h*(r^)r^ (s+ 2 ) 

2 (s+l) (l+K^) (s) 

2(s4-l)sin(vrsMnH-K^J(l+K^) 

|{2s4-3) (l-m) - 

S+l 

M-gh*(-^o)rn 

sli+ryTi+vT 


(s+2 ) cos (tts) 


(m+K )cos ( tts) 


') 


( 82 ) 

(83) 

(84) 

(85) 

( 86 ) 

(87) 

(88) 

(89) 


The integral equations in terms of the unknown density 
functions f(r), g(r)j h(r), and h*(r) are given by the 
conditions that 

00 

T(r)+iN(r) = / [T^j.g(r,r^, 7 r/ 2 )+ iT^gg(r,rQ,7r/2)]drQ , ( 90 ) 
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No(r)=/T^gg(r,ro,7r)dro , (91) 

Nq (i^)=/'r2Q0 O)dro • (92) 

o 


Equations (90) and (9I) are exactly as ( 22 ) and ( 23 ), and 
(92) is due to the presence of the crack on 0 = 0 . N(r) 
and T(r) are the normal and shear stresses on the inter- 
face, and N (r) , N*(r) are the normal stresses on the 
symmetry line in materials 1 and 2 respectively. 

On noting as in Part I that 

f(r^) = g(r^) = 0, c < r^ < 00 

h(r^) = 0 , 0 <r^<a, b<r^<oo 
h*(r^)= 0, 0 < r^ < a*, b* < r^ < 00 


and following the same procedure of Part I leads to the 
three integral equations 

^ / 0(t) - 70 (y)= [Q(y)-iP(y) ] + I 0(t)G(t,y)dt 


+ I e(t)G*(t,y)dt, -1 <y< 1, (93) 


-X 


^ 0 (t) 0 (t)H(t,x)dt 


1+K, 


-1 

1 
-1 


( 9 ^) 


= ~ 2 ^ ^ e(t)J*(t,x)dt-/ 0 (t)l(t,x)dt,-l< x<l , 


1 1 g(t) 1 1 

i t^ + /^|(t)H*(t,x)dt = -2JJJ P*(x)-/ ^(t)j(t,x)dt 


1 
-1 
1+x, 


- f <t> (t)l*(t,x)dt, -1 <x<l , (95) 


-X 


t-x 



^7 


specifying , ^(t) and 4(t). In the following equations, 

if a, a* = 0 then a^,aj is replaced by b,b* respectively. 


H*(t,x) = - 


2tt 


m r 8s^ 12s 1 

l+^£ L(n+s)® in+sfJ 

in 


+ 


r nr -^1 3(1 

_ 1 - 
m 

1 

L ^ + Ki 

, m 1 1 + 

fi - 

m 

T^+S 


* ^ ( It ^2 ) 

I (t,x) = 2 (l+mKj_) (nH-K2)7T 


+ ct 


K, 




-is 


(1 + (5c2t^+ S^) 


c^t^+ s 


_2 


_2 , 


C®t^+ S^ 


J(t,x) = 


7T 


0 


(1+ V) - 

c t 

- J 1 

2 (irH-Kg) 2 ( 1+mx^ ) J ti+s 


3(1+K2) (l+K^) 


c^t^t s^ 


ri+K2 

l+Kg- 

S 

(_ l+IUK^ 


(t]+s)^ 


J*(t,x) = 


3(1+K^) (1+K3_) 


"" 1 L2(^-*-^i) 2(1+^)J n+s 


+ 


u, a* 
^2 o 


" Otl(l+K J tT 


/-iri 


1 + 1 + 


1 + 


Ti " 

m 




(fi+s) 



+ 


K2 

m “ ^1 


_2 2 2 ' 1 
T1 +C y ^ +K^ 


+ cy 


2(!i^-c=Y^)- 


(1- i) (fi^-3c V 


K, 


(•^- + it) 


2(fj^-cy )- 


fi^+c^y^ 


K. 




- K1 


m 


(1- h) (3T)^-cy 

Kg 

+ ~m) 

1 


fi^+c^y^ 
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with ^(tjx)and H(tjX) given in Part I and a* = (b*-a*)/ 2 , 
b* = (b*+a*)/2. 


r a*t+b* , rf a* 7^-0 
loo 

Y j if a*= 0 


a*x+b* 
o o 

b*x 


if--a^ f- 0 
if a*= 0 


Note that the above integral equations could be written down 
immediately by superposition, noting the solutions in [2] 
and Part I. That is, to obtain Equation ( 93 ) from ( 35 ) 
a term must be added to account for the normal and shear 
stress on Q = — due to a crack on 0 = 0 . But whether or 
not the crack is on 0 = 0 or 0 = tt is a matter of nomen- 
clature, and the te]rms corresponding to 0 = 0 may be ob- 
tained from the terms corresponding to 0 = tt by careful 
inspection. To obtain ( 94 ) from (36) the effect of the 
crack on 0=0 must be accounted for, and its effect is 
given in [ 2 ], The third integral equation is given in [ 2 ] 
with the added term given by Part I to account for the 
effect of the interface crack on the symmetry line. 

The solution for </>(t) from ( 93 ) ^ follows directly 
from Part I, and the . are given by 


c = — — V'(t)G(t,y)dt 

L(n-l,-a,-p -iL“i -1 


(-ct,-p) 


+ 


f I (t)G*(t,y)dt J (l-y)~“(l+y)~^Pj^_j^(y)dy 


(96) 


where the term is determined as in Part I, and is given 

by the following: 

if a* ^ 0, a ^ 0 C = 0 

o 

if a* ^ 0,' a = 0 = - ibcosh . ( 7 T , cd ^(t)dt 

O C7T J r \ / 


-1 
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If a* = 0, a ^ 0 




ib*cosh ( , 710 }) 

C 7 T J S V ; 


If a* = 0, a = 0 


C„ = - sj{t)it 

ib*cosh (tto)) 


+ 


C 7 T 


/ I (t)dt . 


-1 


The displacements are then found by integrating (71 )j 
( 72 ), ( 73 ) j from which the constant has been found, 

and the continuity condition of the functions ^(t) and 
4 (t) are specified. Integration of ( 7 I) gives the opening 
displacement on 0= 0, where if a*= 0, a*, is replaced by b*. 


af 1 


u^g(x,o) - -f 4 (t)dt--^/ ^ (t)dt 


( 97 ) 


may now be substituted into (94) and ( 95 ) > and letting 
Q(y) = 0, P(y) = - , one has 

^ /Vt) ^ + / ^(t)k(x,t)dt = P^(x) 

2 Q 

I w(z)l(z,x)p^?^^^dz -/ 4 (t)R*(x,t)dt , 

-1 ' -1 

(98) 


and 


i /' 4(t) ^ + A(t)k*(x,t)dt = ^ PJ(x) 


2 a 1 (a,p) 1 

I w(z)l (z,x)p (z) dz - / ^(t)R(x,t)dt 


a 1-72 -1 


-1 


where 


( 99 ) 


R*(x,t)= J*(t,x)+ - 4 . b*c9 , sh(7rco) 

Pi °° 1 1 (a.P) 

Vi_-y 2 n=l ^ -X 

1 . -P) 

/ G«(t,y)(l-y) “(1+y) %.Jy)dy , 


+ 
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(x,t) = A(z)I*(z,x)dz 

“1 

+ J= 2 L(n-l.-a - ;-p) X 

Vl-y^ n=l 1 r_ry -ft') 

-a,,. 


I G(t,y)(l-y) "^(1+y) ’ 

-1 


k*(t,x) = H*(t,x) + f w(z)l*(z,x)dz 


-1 




+ -= s I,(n_i-_a,-p) J' w(2)I*(z,x)P„(z) dz X 

Vl--v2 n=i '' j ^ 

■^ / /_/^r _Q 


/ G*(t,y) (l-y) (1+y) ^Pj^_i(y)dy , 


and k(t,x) is given in Part I. 

Equations (98) and (99) may be solved in a similar fashon 
to (61) following Part I with the only change being 
rather than having one unknown function to be solved 

for at k discrete points, two unknown functions may be 
assumed F(t^), F*(tj), whose solutions are obtained from 
2 N simultaneous equations that are, in general, not sepa- 


rable. 

For 

a > 0 

let 

fit) 

1 1 

- Fft) with f f(t)dt - 0 

Vl-^ -1 

For 

a = 0 

let 

f{t) 

- ^ F(t) = ^ F(t) 

-vAi-t V i-t^ 

For 

a*> 0 

let 

e(t) 

1 1 

= — F*(t) with / ^(t)dt = 0 

\/l-t^ -1 

For 

a*= 0 

let 

e(t) 

= ^ F*ft) = — - — F*(t) . 

Vl-t Vl-t2 
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Then for example when a > 0 , a* > 0 Equations (97) and 
( 98 ) may be written as 

-Hrk (x 




1 +K, 


ka 


o, - 


^ a^i-y^ -1 


O * V 

/ w(?)l(z,x )P (z) dz , 

2 1 1 


N 


with 

and 
n 


2 ^ F(t ) = 0 , 

k=i ^ ^ 


( 100 ) 


2 1 P(t^)R(x,,t^) + 


k=--i 


r^ 




1+K 


2 a 


(a,^) 


2 / w(t)l*(z,x )P (z) dz , 


t^V 1-^ 


N 


with 


TT 


Z ^ F^(t^) = 0 
k=i 




^k-l)?:" ) 

2 N J ^ 


( 101 ) 


^r = =os(^) 


where r = lj 2 ,...,N-l , tj^ = cos 
and the upper limit in the first and second sum appearing 
in Equations (lOO) and (lOl) has been chosen to be identical 
for convenience. 

Had a been equal to zero and a* been chosen greater 
than zerOj Equations (lOO) and (lOl) would appear as follows, [13] 


N 1 . 

^ 2N + 1 
k-i 


1+K, 


•j— + vrk (x 
. k r 


1 ^ 1 
'+) + N 

J :=i 


i F*(tj)R*(Xj.,t.) 


2 a, 


(a.P) 


2^1 ^ " 


— — / w(z)l(z,x )P (z) dz (102) 

r_ ^ 1 ^ 4 - 

Y 


CL ~ ^ 
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? 2N + 1- 
k=i 3= 


N 1 r 1 - 1 

Z ^F*(t.) ^+7Tk*(x t.) 

HI ^ Lt,-x^ ^ => J 


1+K. 


2a. 


3 ^ 

(a.P) 


2 p, 2 2 


a 


^— / w(t)l*(z,x )P (z) dz , 

2^ _i -^1 


with 


a^s/l-y^ -1 

I J F*(t.) = 0 
j==l " J 

7rr 


(103) 


where r = 1,2,...^N, x^ = cos , tj^ = cos 


r = 1,2,...,N-1, x^ = cos(^) , tj = cos 




and again the upper limit sums on 
assumed equal for convenience. 

The stress intensity factors for 
tions (Sk) and ( 65 ), and for a = 
For a* )> 0 


t, and t . have been 
k 3 

a y 0 are given by Equa- 
0 by Equation (68). 


2[i 

k(a^)= — F*(-l) (104) 

2 

2(x^/a* 

k(b*)=--^ — F*(+l) ( 105 ) 


and for a* = 0 

2M./b^ ^ 

k(b*) = -- ^^^ F*(+l) . (106) 

2 


The strain energy release rate is given by Equation (69). 



, Numerical Solution and Results 

The humefical solution to Equations ( 98 ) and ( 99 ) has 
been obtained for the particular case of aluminum- epoxy 
half -planes for ease in comparing with the results presented 
in Part I. As in Part I, the results presented are for 
plane strain and constant pressure on the crack surfaces 
with no loads at infinity. 

Of particular importance in the investigation is the 
potential for further extension or arrest of a crack origi- 
nating in one half-plane^ as the crack reaches the inter- 
face. In Part I, the instances of a crack reaching the 
interface and either spreading along the bond or crossing 
(without spreading along the interface) was considered. 

It is seen in Figure 11, that extension of the crack into 
the adjacent epoxy half -plane decreases the tendency for 
interface crack growth compared with the results of Figure 4 
and for larger values of half-length c, the strain energy 
release rate is less than for the no-interference interface 
crack. 

Table 3 gives some numerical values of the stress in- 
tensity factors, strain energy release rate, and displace- 
ments for the problem of Figure 11. As the crack extends 
into the epoxy half-plane, k^ initially decreases and k^ 
changes sign as the direction of shear changes. 
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Figure 12 depicts the change in strain energy release 
rate for different locations of a unit length perpendicular 
crack, and, as suggested by Figure 4 of Part I, the geometry 
having the most potential for failure is for the crack to 
intersect the interface from the lower modulus side without 
crossing. 

It seems reasonable to conclude from the results of 
Part I and the present study, that the most critical state 
for continued crack growth occurs when a perpendicular crack 
intersects the interface crack but does not cross. Further, 
for a crack approaching an undamaged interface from the 
higher modulus side, interface damage is much less likely 
than continued extension into the adjacent material. If 
the crack originates in the lower, or equal modulus half- 
plane, the potential for extension, either along the inter- 
face or crossing the interface, exists and depends on the 
relative bond strength and fracture properties of the half- 
planes. The combination of interface growth and extension 
into the higher modulus side appears unlikely. 












^U/Sc* 



Table 3. Numerical Values for Aluminum-Epoxy Half-Planes with a 
b = 1.0, c = 0.1, = 0.4285, cj = 1.0, o = 0.26844. 





































































































AN EXPERIMENTAL INVESTIGATION OF 
BONDED ELASTIC HALF -PLANES WITH AN^. INTERFACE- CRACK 
AND A PERPENDICULAR INTERSECTING CRACK 
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James G. Goree 
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Mechanical Engineering 


James O. Feemster^ Jr. 

Graduate Assistant M. S. Candidate in 
Mechanical Engineering 

In this study, the nature of crack growth near a bonded 
interface between two linearly elastic half-planes is in- 
vestigated experimentally. The accuracy of some specific 
analytical solutions previously developed is demonstrated. 

Material pairs of Plexiglas-Buterate and Plexiglas- 
Plexiglas were used with the interface bond formed by 
using ethylene dichloride. An initial perpendicular crack 
in the Plexiglas and/or an interface crack was developed 
by fatiguing the specimen until the desired configuration 
was reached. The static stress field necessary to initiate 
fracture was monitored with a network of strain gages in 
the vicinity of the crack or cracks and the behavior compared 
with the corresponding analytical solution. 

The stress required to initiate fracture in the per- 
pendicular crack as the near tip approaches the interface 
is shown to be in close agreement with the corresponding 
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analytical solution. The nature of the interface bond and 
the fracture toughness of the Buterate were both such that 
quantitative results could not be obtained as to the effect 
of the intersecting crack on the interface, the interface 
crack, or the perpendicular crack in the Buterate. However, 
some important observations can be made as to the behavior 
of these phenomena. When the perpendicular crack is formed 
in the higher modulus material, the crack tip has a tendency 
to extend into the adjacent half-plane and to leave the 
interface undamaged. In the Plexiglas-Plexiglas material 
pairs, the crack may either cross into the adjacent material 
or extend along the interface, depending upon the relative 
strength of the interface bond as compared to the fracture 
toughness of the two materials. In no instances did the 
crack extend along the interface and into the adjacent 
half-plane simultaneously. 
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INTRODUCTION 

The purpose of this study was to investigate experi- 
mentally the nature of crack growth near a bonded interface 
between two linearly elastic half-planes and to demonstrate 
the accuracy of some specific analytical solutions developed 
for geometries amenable to experimental studies. It was 
assumed that crack propagation in the experimental work 
corresponded to brittle fracture^ according to the classic 
Griffith crack criteria [1], [2].^ 

In the existing analytical solutions appropriate to the 
present work, the specific geometry considered was the case 
of two bonded elastic half-planes containing a finite length 
crack in one material, with the crack being perpendicular 
to the material interface as shown in Figure I 3 , The be- 
havior of the stresses, stress intensity factors, and 
strain energy release rates were investigated, therefore 
indicating the potential for continued crack growth. Solu- 
tions for the crack terminating at the interface or cros- 
sing into the adjacent material have been presented by 
Erdogan, Cook, and Biricikoglu in [3]j [^]» The instance 
of one end of the crack approaching the interface and 
debonding along the interface or debonding and crossing into 
the adjacent plane has been investigated by Goree and 
Venezia in [5]j [6]. 

^ References for this section of the report are listed 
seperately from those of Parts I and II. See page 89 
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SELECTION OF TEST MATERIALS 

It was desirable to select materials with sufficiently 
low fracture toughness so as to allow fracture initiation 
at net stresses well below the elastic limits of the two 
materials. Polymers were chosen as test materials instead 
of metals j because the load levels needed to initiate frac- 
ture would be loWj and polymers could be chemically bonded 
easily. 

After investigating all of the commercially available 
polymers to use for experiments, a material pair of Plexi- 
glas and Buterate was finally chosen. These materials have 
a common solvent^ ethylene dichloride ^ thus allowing them 
to be edge-bonded to produce the desired type of interface. 
Plexiglas was found to be the most suitable of the materials 
available, as it behaves elastically at relatively low 
stesses, and brittle fracture can be initiated easily. It 
was desirable to find another polymer which was dissolved 
by the same solvent as Plexiglas, which had a different 
value of Young’ s modulus, and was not fracture tough. 
Buterate and Polycarbonate were the only readily available 
choices, even though they did not meet all of these re- 
quirements. Although they were the only other polymers 
available that use the same solvent as Plexiglas, both 
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materials were extremely fracture tough. Buterate was chosen 
instead of Polycarbonate ^ because a larger modulus ratio 
could be achieved with the Plexiglas-Buterate material pair. 
Because of the resistance of Buterate to fracture^ it was 
used as the right half-plane in all of the experiments in- 
volving Plexiglas-Buterate material pairs. Providing that 
fracture was initiated in the Plexiglas, the fracture tough- 
ness of Buterate did not affect the outcome of the experi- 
mental work. However, fracture along the bonded interface 
was influenced by the toughness of the Buterate. 

Another consideration in the choice of materials was 
the thickness of each test specimen. It was decided to 
use 1/8 in. (.3175 cm) thick sheets of Plexiglas and 
Buterate, as this thickness would be small enough to ap- 
proximate the desired state of plane stress, while being 
thick enough to provide sufficient surface area to give a 
satisfactory interface bond between the two materials. 
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EXPERIMENTAL PROCEDURE 


All experimental work was conducted on an Instron Model 
109 dynamic testing machine as shown in Figure 14 „ With the 
aid of a network of strain gages placed on each test sped-- 
men, the stresses exerted by the machine could be monitored. 
Of considerable importance was the ability to obtain a stress 
field near the cracks which closely approximated the uniform 
stress field assumed in the analytical solution. Therefore ^ 
it was desirable to make each test specimen as large as 
possible within the limits imposed by the size of the test- 
ing machine, in order to minimize the effects due to the 
machine grips. Each material pair tested was 12 in. ( 30.48 
cm) in height and I5 in. (38. 10 cm) in width. 

In order to calculate the stress in each test specimen, 
it was necessary to determine experimentally the value of 
Young’s modulus E and Poisson’s ratio v for each material. 


The following values were found; 


E = 
1 


'Plexiglas 


= 5.07 X 10 ^ psi 




■'’Plexiglas = , 


LL„, . , = 1.856 X 10 ’= psi ( 1.280 X 10 ® 

’^Plexiglas \ 


'Buterate 


i (1. 


= 2.58 X 10 ^ psi (1.780 X 10 ® ^) , and 


V = V 
2 


^ ^ = 0.37, LLo ^ --::9.4l6xio^psi(6.496 xio® 

Buterate ^Buterate ^ \ niy 



65 


In the analytical solutions [5]^ [ 6 ]^ the material 
half-planes jvere assumed to be loaded with uniform stresses 
o j and a as shown in Figure 13, with the strains 

O 1 2 

being related in such a manner as to give constant strains 
in both the X and Y directions at points remote from 
the cracks. The following relations between the stresses 
must then hold: 



for generalized plane stress, and 



for plane strain. 

In order to achieve uniform strain in the Y-direction 
across the specimen shown in Figure 13, the resultant load 
must act at a distance d from the material interface, as 
shown in Figure 15 . This distance d for a specified 
material pair can be computed by first finding the rela- 
tionship between the stresses and , and then by 

summing moments about the material interface. By refer- 
ring to Figure 15 for a specimen of thickness t, resultant 
loads are 

Ri = tw^ , and (3) 

R = cr tw 
2 22 


( 4 ) 
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Summing moments about the interface of the specimen^ 


Rd = 




R w 
2 2 


where 


( 5 ) 


R = R + R 

X 2 


( 6 ) 


Substituting Equation (6) into Equation (5) gives the fol- 
lowing relation for d : 


d = 


2 (a w + a w ) * 
11 22 


( 7 ) 


As an example j the distance d for a material pair of 
Plexiglas-Buterate in the state of plane stress with 


w = 


w = 


6 in. (15.24 cm) j 
6 in. (15.24 cm) j 


a = 1.0^ and 


o_ = 


0.5 . 


will be calculated. 

Using the experimentally determined values of E and v for 

Plexiglas and Buterate, a is given by Equation (l) as 

2 

O = 0.6007 o , 

2 1 

and Equation (7) then yields 


d = 


UlWi^ - ^2^2^ 

2(cJ^w^+ o^w^) 


= 0.748 in. (1.899 cm) 
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liT the resultant R of the load is placed this distance from 
the interface in the higher modulus side of the material ~ 
pair^ uniform strain across the specimen in the Y-direction 
will result. 

In the adaptation of the analytical problem to experi- 
mental workj it was found that this distance d could be 
approximated as zero^ due to the manner in which each 
specimen was constrained during loading, and uniform strain 
parallel to the interface would still be obtained. Figure 16 
shows the method of gripping each specimen for testing. The 
grips were fabricated from eight steel bars 3/8 in. (. 952 cm) 

X 1 in. (2.5^0 cm) X 12 in. (30.480 cm), which were bolted 
tightly on each test specimen, two on opposite sides of 
each of the four edges of the specimen. Underneath each 
steel grip was placed a strip of Plexiglas 1/8 in. (.317cm) 

X 1 1/2 in. (3,810 cm) X 12 in. (30.480 cm), in order to 
reinforce each test specimen and to distribute the load 
more evenly. 

By considering the actual method of loading each 
specimen, the value of d, as shown in Figure 16^ will 
again be calculated. Using the same Plexiglas-Buterate 
material pair in plane stress with 


w^ = 6 in. (15.24 cm), 
w^ = 6 in. ( 15.24 cm), 


a 


0 


1.0 


9 


and 
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Oq = 0.5 , 

the value of d calculated will be shown to be negligible.^ 
For the purpose of computing this distance, the resultant R 
is shown in Figure 16 to be placed a distance d from the 
interface in the higher modulus material. It should be 
noted here that the stress field near the grips is not uni- 
form. However, in order to obtain an approximate value for 
the distance d , which accounts for the presence of the 
lateral steel grips, the strains parallel to the interface 
will be taken as equal in all three materials and given by 
the average strain in either the Plexiglas or the Buterate, 
as computed from the uniform stress fields or . 


Therefore, 


^Plexiglas ^Buterate ^Steel 


The relationship between the stresses and in the 

Plexiglas and Buterate was found previously to be 

a = 0.6007 cf 

2 1 

The value of the three strains can be found as follows: 


^Buterate Eg ^^2 '^s^o^ ^Plexiglas ^Steel 


where 


= 2.58 X 10 ^ psi (1.780 X 10 ® ^ , 


and 


= 0.37 . 



Then^ 
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^Buterate 2.682 X 10 in/in (2.682 X 10 . 

The approximate tensile stress in each of the four lateral 
steel grips can now be found; 


= e 


Steel Steel 


= 8o. 46 0 


where 

e„. , = 2.682 X 10 ^ 0 in/in (2. 682 X 10 — ),and 

steel 2 2 cm 

Esteel = 30 X 10® psi ^2.070 X 10® ^ 

The resultant R in Figure l6 is now found to be 

R = 6 0^t + 6 0^t + ^ o'g (l) (3/8). 

Summing moments about the interface of the specimen j 


Rd = 6at-6a t . 

1 2 

Then^ 

6ot - 6o^t 

d = -T — ■ , ^ , — ■p — = 0.004 in. (O.OlO cm) . 
oa t + 6a t+ 1.50 

12 s 

Since this value of offset is then effectively zero 
within the accuracy of the tests, the resultant load under 
actual test conditions was placed along the centerline of 
each specimen. During testing, the condition of uniform 
strain across the specimen parallel to the interface was 
verified with the use of strain gages. 

For each specimen tested, it was desired to establish 
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a specified ratio between the stress cr and the lateral 

1 

stress . This objective was achieved with the use of 

eight alien head screws threaded into the side grips as 

shown in Figure l6 . By tightening these screws against 

the top and bottom grips ^ the side grips could be forced 

outwardj thus creating a lateral stress on the specimen. 

It should be noted that high values of a , such as 

= 1.5cr^ j could be easily achieved with this method. 

In order to establish the desired values of a and a , 

1 o 

the strain in the X and Y directions was carefully moni- 
tored by a network of strain gages, as is shown in Figure 16. 
Strain gages having 350 ohm resistance were used, because 
the problem of gage "drifting" associated with heating 
could be eliminated. Gages having 120 ohm resistance, and 
therefore larger current, had a tendency to store heat, 
because the polymers that were used as test materials were 
not good thermal conductors. In the initial tests, twenty 
strain gages were placed on each test specimen in order to 
investigate the stress state in the vicinity of the cracks. 
The information gathered in these initial tests suggested 
that the gage placement as shown in Figure l6 was both 
economical and practical, and that the stress state moni- 
tored in the vicinity of the cracks was uniform to within 
four percent. 

In order to achieve satisfactory correlation between 
the analytical solutions and the experimental results, it 
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was necessary to form the cracks by fatigue to a desired 
length in each test specimen. With the use of the Instron; 
Model 109 dynamic testing machine^ the mean load, amplitude, 
and frequency of loading could be specified, thus enabling 
the operator to fatigue the test crack to any desired 
length. 

Before the tests on the bonded material pairs were 
conducted, the critical stress intensity factor for 

Plexiglas was determined. In order to find this value , 
five test specimens were made and tested, as shown in 
Figure I 7 . In each case, a 5/32 in. (0.397 cm) hole was 
drilled in the specimen, and a notch was made with a 
surgical scalpel on each side of the hole perpendicular 
to the direction of loading. A crack was first fatigued 
to a desired length, and then loaded statically to de- 
termine the critical stress o at fracture. Strain 

1 

gages were used to accurately monitor this stress. The 
value of was then deteimiined from the following rela- 

tionship: 

K = a , 

^ X 

where 

a = stress at which the crack propagated, and 

t = half-length of the initial crack . 

The value of for the Plexiglas used in the experimental 

work was found to be 

Kc = 530 psi Vin. (5828 ^ Vein ) o 
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After determining the value of for Plexiglas, tests 
were then conducted on the Plexiglas-Buterate^ and Plexiglas- 
Plexiglas material pairs. First, the materials were edge- 
bonded to form an interface, and reinforced around the four 
edges with eight Plexiglas strips 1/8 in. (0.317 cm) X 1 1/2 
in. (3.810 cm) X 12 in. (30.480 cm). In each case a hole 
1/4 in. (0.635 cm) in diameter was drilled in the left half- 
plane a specified distance from the interface. A scalpel 
was again used to make a notch in one side of the hole 
nearest the bonded interface, as shown in Figure I6. The 
notch was made on one side only in order to keep the crack 
from propagating to the left under the fatigue load. If an 
interface crack in the test specimen was desired, a 1/32 in. 
(0.079 cm) wide slot was milled a desired distance along 
the interface, and a scalpel was used to put a notch at 
each end of the slot. At this time, the network of strain 
gages was placed on the specimen. The specimen was then 
secured in the testing machine, and turned perpendicular to 
the usual axis of loading, and an interface crack was fa- 
tigued. In all* of the tests performed, a fully satisfac- 
tory interface crack was never achieved, as the solvent 
used in bonding the materials changed the material proper- 
ties slightly, thus making them more fracture tough. How- 
ever, a more satisfactory interface crack was obtained in 
the Plexiglas-Plexiglas material pairs than in the Plexi- 
glas-Buterate combination. 
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The test specimen was next placed in its normal loading 
position, as is shown in Figure 16 , and was fatigued under 
a specified load and frequency until the end of the crack 
perpendicular to the interface had propagated to within a 
desired distance of the interface. It was then loaded 
statically under a specified ratio of a^/a^ , which was 
monitored by the network of strain gages previously de- 
scribed, to determine the stress at which the end of the 
crack started to grow toward the interface. The effect 
of the crack upon reaching the interface or the interface 
crack was then observed. It should be noted that the 
presence of the drilled hole at the left end of the crack 
did affect the outcome of the experimental results for 
short cracks, as is shown in the following section. 
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EXPERIMENTAL RESULTS AM) CONCLUSIONS 

In Figure l8, a comparison of the analytical and ex- 
perimental results in Plexiglas-Buterate material pairs is 
shown by curves 1 and 2 respectively. The critical stress, 
which is the stress required to initiate fracture, is plot- 
ted versus the distance of the crack tip from the inter- 
face. According to the results of [3] through [6], the 
critical stress approaches zero as the crack tip nears the 
interface, and approaches infinity as the crack length be- 
comes zero. Curve 3 of Figure 18 gives the critical stress 
required for a comparable length crack in a full-plane of 
Plexiglas. The critical stress for the full-plane problem 
approaches infinity as -t approaches zero, as is shown 
from the following relationship: 

a = K / 
c c 

where a = critical stress 

c 

= critical stress intensity factor, and 
-t = half-length of the crack . 

The experimental results agree closely with the 
analytical results in the region a = 0.05 in. ( 0.127 cm) 
to a = 0.4 in. (1.016 cm). However, in the region where 
the crack tip is in the vicinity of the drilled hole, the 
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stress concentration associated with the presence of the 
hole reduces the critical stress , as shown in curve 2 
of Figure l 8 . The problem of determining the stresses and 
stress intensity factors for symmetric cracks emanating 
from a circular hole in a full plane was considered by 
Newman in [7]. The results given in TABLE 1 of [7] 
reproduced in Table 4 of the present study. Note the 
column corresponding to in TABLE I of [7] should read 

A = 0 rather than A = 1. Similarly^ the column for F^ 
should read A = 1, These changes have been made in Table 4 
of this investigation. The present one-sided crack seems 
to be influenced by the presence of the hole more strongly 
than the symmetric crack does. However, the trend is the 
same. 

In Figure 18, curves 1 and 3 agree within 5 percent 
of each other until the crack tip is within approximately 
0.2 in. ( 0.508 cm) of the material , interface. The crack 
tip is then affected by the presence of the Buterate in 
the right half-plane, as is shown in curve 1 , and reduces 
the critical stress required for fracture. 

Because of the inability to produce a brittle in- 
terface, no quantitative results were obtained for the 
behavior of an interface crack or the effect of the per- 
pendicular crack on initiating fracture along the bond. 
However, some important observations were made as to the 
influence of the modulus ratio of the materials in both 
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Plexiglas-Buterate and Plexiglas-Plexiglas material pairs. 

In the Plexiglas-Buterate combination, an interface 
crack could not be formed as a result of a perpendicular 
crack intersecting the interface. According to the analy- 
tical solution [3]j a material pair with the intersecting 
crack in the more rigid material actually has compressive 
stresses induced on the interface in the vicinity of the 
crack tip. When the perpendicular crack tip reached the 
interface, arrest occurred in all experiments, regardless 
of whether the test specimen was under a static load or a 
fatigue load at the time of intersection. By continued 
fatigue loading, the crack could always be forced to cross 
into the Buterate, as shown in Figures I 9 and 20 , and af- 
ter a crack growth of approximately 0.1 in. ( 0 . 25^0 cm), ' 
brittle fracture would occur. As discussed in [3]j the 
order of the stress singularity for the cleavage stress of 
the intersecting crack is larger than one-half, providing 
that the crack originates in the higher modulus material, 
and it is assumed that this large opening stress was suf- 
ficient to develop the fatigue crack in the Buterate. It 
is of considerable interest that a fatigue crack could be 
produced, and that brittle fracture would occur in the 
Buterate in this manner. In a full plane of Buterate with 
an initial configuration as shown in Figure I 7 , fatigue 
loading in the identical load and frequency range would 
not generate a crack. After fatiguing for 200,000 cycles 



at a load of 1400 + 200 lbs. and a frequency of 15 hz. ^ 
the crack tips in the full-plane of Buterate became highly 
distorted^ and multiple irregular cracks were formed^ as 
shown in Figure 21 , none of which lead to fracture. With 
a large lateral load, the interface in Plexiglas-Buterate 
material pairs could be made to fail in tension, as it was 
weaker than either of the two materials. In one experi- 
mental run, a crack originated on the interface due to a 
flaw in the bond approximately two inches from the inter- 
secting crack tip. As the lateral load was increased, 
this crack behaved similarly to a brittle fracture, and 
it propagated along the interface until the intersecting 
crack tip was reached, at which time it stopped, leaving 
the bottom half of the interface still intact. 

Unlike the Plexiglas-Buterate pair, a crack could be 
initiated along the interface in the vicinity of the in- 
tersecting crack tip in the Plexiglas-Plexiglas material 
combination, depending upon the lateral stress, and on 
the intersecting crack length. As in the Plexiglas- 
Buterate material pair, the interface always stopped the 
intersecting crack. Both the lateral stress, , and 
the intersecting crack length were significant in pro- 
ducing interface failure. However, the nature of the in- 
terface bond was such that consistent results could not 
be obtained. With large values of the lateral stress, 
such as a = 1.0 a 




the interface would first arrest 
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the intersecting crack tip, but under continued fatigue 
loading the interface would fail, providing that the in- 
tersecting crack length was 0.8 in. (2.0320 cm) or less. 

For intersecting cracks longer than this value, the crack 
tip would cross the interface, resulting in fracture of 
the test specimen. 

In conclusion, for material pairs having |x /\i <( 1 , 

compressive stresses exist on the interface in the vicinity 
of the intersecting crack. In this case, the strength of 
the interface bond is relatively unimportant, because 
failure is more likely to occur in the form of the inter- 
secting crack crossing the interface into material 2. How- 
ever, according to the analytical solution [3]> when 

|i /[x )> 1 , the interface is subjected to tensile stresses, 

giving rise to both interface fracture or fracture through 
the interface into material 2 as possible modes of failure. 
The type of failure which occurs is then dependent on the 
strength of the interface bond as compared to the strengths 
of the two materials being used. In none of the experi- 
ments involving Plexiglas-Plexiglas half-planes did simul- 
taneous crack extension occur along the interface and 
across the interface. This is in agreement with the 
analytical results discussed in [6]. 
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Figure 13. Bonded Elastic Half-Planes with Perpendicular 
Cracks . 
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Figure 14. Instron Model IO 9 Dynamic Testing Machine. 
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Figure 15. Theoretical Resultant Load Placement for a 
Plexiglas-Buterate Material Pair. 










o 


o 


o 


o 



Figure 17. Configuration of Plexiglas or Buterate Full 
Plane Prior to Fatigue Loading. 
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DISTANCE a OF CRACK TIP FROM INTERFACE (in.) 


Figure 18. Comparison of Analytical and Experimental 
Values of Critical Stress for Plexiglas- 
Buterate Material Pairs, b = 0.8 in. 

( 2.0320 cm), d=0,25 in, (0.6350 cm). Plane 
Stress . 
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Figure I9. Perpendicular Intersecting Crack in a Plexi- 
glas-Buterate Material Pair Crossing into the 
Buterate Half-Plane Following Continued Fatigue 



Figure 20. Magnified View of Crack Crossing into the 

Buterate Following Continued Fatigue Loading 
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Figure 21 


Irregular Crack Formation in a Full-Plane of 
Buterate Following Continued Fatigue Loading 


Cracks Emanating from a Circular Hole 
in an Infinite Plate Subjected to 
Biaxial Stress from [7]. 



F(A= -1) 


F^(A = 0) 


F^(A = 1) 


1.0135 

1.0775 


1.2208 


1.2350 
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. 6082 
.7104 
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.8400 

.9322 

.9851 

1.0168 

1.0358 

1.0536 

1.0582 

1.0571 
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1.0409 

1.0336 

1.0252 

1.0161 

1.0077 
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APPENDIX A 


ON THE TRANSFORM SOLUTION 


1 . The Mellin Transform and its use for the Plane Problem 
of Linear Elasticity 

Following [1] and [14], the plane problem of linear 
elasticity may be formulated as follows: 

Trg + iTgg = a(r, 0 )= ~ ^ ^ ^ 1 * 1 ) 


rr - ^ r Sr ^ 


(A 1.2) 


2|xl 




= v(r e)= - + -4 ^ - i 

Sry ^0 r SrS 0 


+ (l-v)fe + 2 ir If + r^l^ + ir^ 


2 J- 


Sr^J 

(A 1.3) 

where [i is the shear modulus, V is the Poisson’s ratio, 
and 


V*x = 0 , = 0, ^(r|f) = V-x=0 + i|f + 4 


^ S0^* 


(A 1.4) 

The Mellin transform of a function f(r), defined and suit- 
ably regular on ( 0 ,«> ), and its inverse are defined by 

F(s) = M[f(r)] = /“f(r)r®“^dr , (A I. 5 ) 


1 rC+i“>_ , \ -S 


^ . ^(s)r ds , (A 1.6) 


c-100 


where c is such that r f(t) is absolutely integrable 




A-2 


on (0, “) . The transform of derivatives is given by 


j-“ j,s-i+n ^ Ils±^ F(s), (A I. 7 ) 


provided 


o dr 


^ ^ -* 0 as r -*■ (O, , m= 1, . . . ,n . 

dr^ 


(A 1.8) 

The regularity conditions (A 1.8) provide the criterion for 
the choice of c . 

Taking the Mellin transform and solving (A 1.4) from 
(a 1 . 1 ) -(a 1 . 3 ) we obtain 

r i(s+2)e _ -i(s+ 2 )e-| 

M[r^cr(r,e)] = 2(s,e)=2i(s+l)Ase-^'=’'^+B(s+l)e -Be 

(A 1 . 9 ) 


M[r^T (r,e)]=-s(s+l) (Ae^®^+Ae"^®® ) - (s+l) (s+4) X 


r ^ 

Be 


i(s+2)e _ -i 


M[ ^^v(r ^_ e l]^V(s, 0 )= - Ase^^^+B(s+l)e +KBe 

M' ^ 


i(s+2)en 

+ Be J . (A 1.10) 

i ( s+ 2 ) 0 _ -i ( s+ 2 ) e 1 

Be J , 


(A 1.11) 


where k = 3 ~ 4v for plane strain k = (3~v)/(l+v) for 
generalized plane stress , and the "integration constants" 

A and B are complex and are determined from the boundary 
conditions specified along the wedge boundaries, 0 = con- 
stant. After A and B have been determined, it is necessary 
to invert (A I. 9 ) through (A l.ll) to determine the stresses 
as a function of r and 0. The next section deals with 
the evaluation of particular integrals applicable to the 
problem in question. 
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2 . Some Particular Inverse Mellin Transforms 

Given beloWj for reference, is a listing of some of 
the Mellin inverses used in inverting (A I. 9 ) through 
(a 1 . 11 ). The details of a few of the integrals will be 
included for completeness. 


1 

27ri 

1 

27Ti 

1 

2Tri 

■1 

27Ti 


1 

27Ti 


1 

27ri 

. 1 
27Ti 


ctioo „ -s-2 

; (^) ds = ) 

C-lOo O 


c+i“ 

1 

't' 

-s -2 

ds = 


2 

r 

0 

-ioo 

sin (tts /2 ) 


( 


C+ia> 

1 


-s -2 

ds = 


2 r r 
0 

—ioo 

cos (773/2) 


77 (r 

^+r^) 

C+i°® 

s 

(^) 

0 

-s-2 

ds = 

i 


0 

-loo 

sin (773/2) 

77 

(rj+r^)^ 

c+i°° 

s 

0 

-s -2 

ds = 

2 

r 

3 r'^+r^r^ 
0 0 

-ioo 

cos (t 7 s/ 2 ) 

77 

r 

0 

_ (rj+r^)\ 

c+i«> 

sin( 77 s/ 2 ) 


-s -2 

ds = 


2 rr T 

0 1 

-ioo 

cos (t 7 s/ 2 ) 


77 

(r^-r^) 
' 0 ' 

c+i«> 

cos (t 7 s/ 2 ) 

(A 

0 

-s -2 

ds = 


ro 
0 to 

-ioo 

sin( 77 s/ 2 ) 

77 ( 



(A 2.1) 
(A 2.2) 
(A 2.3) 
(A 2A) 

(A 2.5) 
(A 2.6) 

(A 2 . 7 ) 


Consider now the evaluation of (A 2.2) and (A 2.4) 


= f 

2 tti •' 


c+100 


c-100 


-s -2 

(— ) ds , r ^ 


Let I 


sin (tts/ 2 ) 



applying the residue theorem and closing the contour to the 
left as shown, noting the zeros of sin (tts/ 2) are at s=+ 2n, 



^ 00 2n-2 ^ 00 nvfl 2m o °o n 2n 

I = I 2 (-l)’^(f-) = f S (-1) i~) = - f 2 (-1) (— ) 

^ n=l o m=0 o n=0 o 


2 

7T 




2r! 


TT(r^+r^) 


as given by (A 2.2). 


For r y 3 the contour is closed to the right and gives 


the same result. 


For (A 2.4) let I = ^ . 


Otioo 


-s-2 


(”) ds, r < r 


27Ti . . / /„\ 'r ' 

C-loo sin{,7Ts/2) o 

again applying the residue theorem and closing the contour 


to the left. 



g" 


I = |£ (-gn) (-l)“(^) = - ^ S 2(mfl)(-l)"'^(f-) 

n= 1 o n= 0 o 
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2m 


^ m=0 o 


noting that 


“ n r- 2n 1 

2 i-ifif-) =— 7 ^ ^ 

n=0 ""o l+(^) 

r\ 


and differentiating both sides with respect to (■“) gives 

2n-l 


S (-I)"" 2n(^) 


2 if-) 

o 


n=l 


1 +(#-) 
o 


rH-1 , j- ,2n+l 


S2(n+1)(-1) (-^) 


2 (^) 


n=0 


- — 2 (n+1) (-1) (— ) 

o n=0 o 


2 if-) 

o 

■ r y 2-, 2 

o 

- ^ 




2n 


n=0 




and I = i 




as given by (A 2.4) 


For r y the contour is closed to the right and 


again 


gives the same result. 
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A similar approach gives 


1 

27ri 


C+ioo 


/ 


1 

sin (tts ) 


-s-2 

(— ) ds 

\r / 



1 

27ri 


c+i“ 


r 

s 

J 

c-i°o 

sin (tts ) 


-s-2 

(— ) ds 

'r ' 


2r r r 
o o 

7r(r^+r)2 


(A2.8) 


(A2.9) 


c+ioo - . 

1 r COS (tts ) 

27 Ti ■' sin (tts ) 

C-loo ' ' 




(A2.10) 



APPENDIX B 


ON THE METHOD AND ACCURACY OF THE NUMERICAL SOLUTION 


1 . The Method of Numerical Integration 

The numerical integration was performed using four 
basic integration Gaussian- type formula [12, p. 887 - 89 O] 


f ^ 


dx 


X, = cosi|^ 


7T 

w . = — 
1 n 


(B 1.1) 


n 


n 


- l T ■“ O -i ^ 27T (l+Xi ) 

' *i= gn+1 

(B 1.2) 

/ f(x) Vl-x2 t S^w^f(x^), x^ = cos(^), sin^(x^) 

(B 1.3) 

^ f(x)dx s; 2 Wj_f(x^),x^= ^ 

i=l 


where 


y . = X 


th 


zero of P (x), w. 

n ' ^ 1 


(l-yp[p;(Xi)] 


(B 1.4) 


and P^ is the n^^ order Legendre Polynomial. 

Wherever possible advantage was taken of the fact that the 
discrete values of F(tj^), given by ( 63 ) of ( 67 ), were at 
the collocation points of (B l.l) or (b 1.2) respectively 
if n was chosen equal to N. However, in the integrals 
specifying the displacements it was necessary to use integra- 
tion formula whose limits were not, in general, -1 to 1 . 
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For these integrals the discrete values of were 

fitted by least squares analysis to a continuous function 
in the form of a order polynomial and an integration 

formula of the form (B 1.4) was utilized. 

2. The Accuracy of the Numerical Solution 

The accuracy of the computer solution was determined 

in part by to what degree it could reproduce known solutions 

or how well it reduced to special cases whose solutions were 

known exactly. For instance, the special case of a unit 

crack in a full plane could be modeled one of three ways: 

either by taking c=0,V =V ,E =E ,a=l,b=2;or 

1212 

c = 1/2 ,a=0,b=0,v =v ,E =E ;or lastly, it 

12 12 

could be approximated by taking c = 0 , a large, and b = a +1. 
The numerical solution was correct to five significant fig- 
ures for all special cases and was exact to the number of 
places carried, (usually four), for the known solutions of 
a crack approaching a bonded interface given by [2], While 
these checks provide a measure of the accuracy, they are 
also necessary for the correctness of the computer solution. 

The terms that necessarily go out in checking the 
special cases are the additive terms that provide for the 
new solution. To make sure that these terms were programed 
correctly, the computer solution of the integral equations 
was approached two ways, and the programs, independently 
written, were checked until agreement was reached. 
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The method that was not chosen for inclusion in the text was 
an iterative solution obtained as follows: As a first 

approximation it was assumed that , given by Equation 

(5^)j was the only non zero . Knowing the ^ the 

right side of Equation (36) was completely defined and a 
set of could be determined. Based on this set of 

F*s, a new approximation to the could be found from 

Equation (53). The right side of Equation (36) was again 
modified and a new set of F' s determined. The procedure was 
continued until no significant change was determined in 
either the F's or the . The accuracy of the complete 
solution was then determined by noting to what extent these 
two solutions agreed. For example, if M = 14, N = 14, 
a= 0, b= 1, c= 1, the iterative solution converged 
after four iterations, giving k(b) = .9513^ 1.117X10", 

whereas, the computer program, discused in Appendix C, gives 
k(b) = .9511. = 1.116 X 10"® . 



APPENDIX C 


COMPUTER PROGRAMS 


1 Introduction to the Computer Programs 

Two computer programs were written to obtain numerical 
results to the solutions presented in Chapters I. and II . 
These will be referred to as PROGRAM 100 and PROGRAM 200 
respectively. Adequate comments have been provided for the 
reader familiar with Chapters I . and II to follow the logic 
of the programs. Wherever possible ^ these comments refer 
directly to equations in the text. Given below are two 
definitions that may prove helpful in the reading of the 
computer programs. The first is a fundamental integral 


containing which was indexed to limit the use of the 

recursion formula used to generate the Jacobi polynomials. 

1 (a,P) 


PLSTAR(N,KX) ^ / w(t)P (t) I(t,x, )dt 

-II K 


(C 1.2) 


The second integral contains P^ ^ (y) . It is related to 
the C^ and is given as part of Equation ( 53 ). Its defini- 
tion follows from 


KQUAD XI p ^ 4- ^ 1 

S F(TJ)PMSTAR(N,J) ^ J f -lUl s(t,y) X 

j=l -1-1 V 1-t^ ^ 

(-Op-P) 

(c 1.2) 


where F(t) is given by Equation (63) and KQUAD is a 
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numerical constant the choice of which is discussed in the 
section that follows. 

2 . On the Choice of Constants Related to the Numerical 
Integration 

. Certain constants have been fixed in programs 100 and 
200. The choice of these constants was governed by how a 
variation in the parameter effected the most sensitive 
output variable in trial runs using Alviminum-Epoxy material 
pairs. The constants along with their definitions are 
listed below for reference. 

KQUAD - Upper limit in numerical integration formula 
used in integrating Equations (53) and (6o). 

Fixed in program at value of 200. 

NTERMS - Number of terms in polynomial curve fit used 
in numerical integration formula (B 1.4). 
in program at value of 10, 

NODE - Number of nodal points used in evaluation of 
displacement integrals that were complete 
and could be evaluated using integration 
formula (B l.l) or (b 1.2). Fixed in 
program at value of 100. 
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3 PROGRAM 100 

Purpose 

The solution is given for two bonded isotropic linearly 
elastic half-planes of different elastic properties having 
a crack along the interface as well as a perpendicular 
crack in one of the half-planes which may intersect the in- 
terface crack. A constant pressure is assumed on the crack 
surfaces j and no loads at infinity. For convenience , the 
inplane crack will be assumed to be in material one for all 
cases. The geometry is shown below. 
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DESCRIPTION OF PARAMETERS 


NSETS Number of data sets to be read, a data set 

consists of cards No. 2 to No.5j as described 
in INPUT. 


A 

B 

C 

NUl 

El 

NU2 

E2 

SIGMAO 

SIGMAl 

NHALF 

N 

M 


STRAIN 


INPUT 

DATA CARD 


Distance of near crack tip from interface. 

Distance of far crack tip from interface. 

Half length of interface crack 

Poisson’s ratio, material one. 

Young's modulus, material one. 

Poisson's ratio, material two. 

Young's modulus, material two. 

Opening pressure on interface crack. 

Opening pressure on inplane crack. 

Number of half planes present. 

Number of terms taken in Equation ( 63 ) 
and ( 67 ) (11 < N < 28) . 

Number of C^ generated by Equation (53) 
(1<M<50). 

If strain set equal to 1, solution is for 
plane strain. If strain set equal to 2, 
solution is for generalized plane stress. 


PARAMETERS 


FORMAT 


No.l 

NSETS 

12 

No. 2 

A, B, C 

3(IPE14.7) 

No. 3 

El, NUl, E2, NU2 

4 (IPE14.7) 

No. 4 

SIGMAO, SIGMAl 

2(IPE14.7) 

No. 5 

N, M, STRAIN, NHALF 

4(12) 


OUTPUT 


Numerical results are obtained for the stress intensity 
factors j strain energy release rate, stresses and displace- 
ments . 

Comments 

In test programs using Aluminum-Epoxy and Epoxy- 
Aluminum ma,terial pairs, the output was examined using 
successively larger values of M and N to determine when the 
numerical solution was sufficiently conditioned to give 
less than 1 percent variation in the most slowly converging 
output parameter. 

It was found that the strain energy release rate con- 
verged most slowly in all cases and that K(b) converged 
most rapidly in all cases. Values that gave less than 1 
percent variation in output atA= 0, B= 1, C= 1, were 
M = N = l4, If A = 0, B = 1, C = 0.1, it was necessary to 
take M = 30 and N = 20. The value of N = 26 was found to 
be large enough in all cases. If C is decreased below 0.1, 
it was necessary to take larger values of M to obtain 
accuracies of one percent in the strain energy release rate. 
PROGRAM, 100 will solve the following problems illustrated 
below. Note A and C may not equal zero at the same time. 
Time required on the Clemson University IBM 370/158 computer 
was about two minuets for N=M=l4 and fifteen minuets for 

n=26, m=48 . 
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A^O, B^O, C=0, NHALF=2 


A^O, B^O, C^O, NHALF=2 


A=0, B^O, Ci^O, NHALF=2 



A = 0, B^O, C= 0, NHALF=1, 
E2, NU2 any nonzero con- 
stant 



A^O, B^Q, C=0, NHALF=1, 
E2, NU2 any nonzero con- 
stant 
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4 . PROGRAM 200 

Purpose 

The solution is given for two bonded isotropic elastic 
half-planes of different elastic properties having a crack 
along the interface as well as a perpendicular crack in 
both of the half-planes, either of which may intersect the 
interface crack. A constant pressure is assigned on the 
crack surfaces and no loads at infinity. The geometry of 
the problem is shown below. 
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DESCRIPTION OF PAJIAMETERS 


NSETS 

A 

B 

C 

AS TAR 

BSTAR 

NUl 

El 

NU2 

E2 

SIGMAO 

SIGMAl 

SIGMA2 

N 


Niunber of data sets to be read. A data 
set consists of cards No. 2 to No.5> as 
described in INPUT. 

Distance of near crack tip from interface 
in material one. 

Distance of far crack tip from the inter- 
face in material one. 

Half length of interface crack. 

Distance of near crack tip from inter- 
face in material two. 

Distance of far crack tip from interface 
in material two. 

Poisson’s ratio, material one. 

Young’s modulus, material one, 

Poisson’s ratio, material two. 

Young's modulus, material two. 

Opening pressure on interface crack. 

Opening pressure on inplane crack in 
material one. 

Opening pressure on inplane crack in 
material two. 

Number of terms taken in Equations (lOO) 
to (103) (11 < N < 28). 


M 


Number of generated by Equation ( 96 ) 
(1 < M < 50). 


STRAIN = 1 implies plane strain, 

STRAIN = 2 implies generalized plane stress. 


STRAIN 



INPUT 


DATA CARD 

PARAMETERS 

FORMAT 

No.l 

NSETS 


12 

No. 2 

A, B, C, 

ASTAR, BSTAR 

5(1PE14.7) 

No. 3 

El, NUl, 

E2, NU2 

4(1PE14.7) 

No. 4 

SIGMAO, 

SIGMAI 

3(lPEl4.7) 

No. 5 

N, M, STRAIN 

3(12) 


OUTPUT 

Numerical results are obtained for the stress intensity 
factors, strain energy rate, stresses and displacements. 


Comments 

In test programs using Aluminum-Epoxy and Epoxy- 
Aluminum material pairs, the output was examined using suc- 
cessively larger values of M and N to determine when the 
numerical solution was sufficiently conditioned to give 
less than 1 percent variation in the most slowly converging 
output parameter. It was found that the strain energy 
release rate converged most slowly in all cases. K(b) and 
K(BSTAR) were the most rapidly convergent output parameters 
in all cases. Values that gave less than one percent 
variation in output atA=0, B= 1, C= 1, ASTAR=0, BSTAR=1, 

were M = 1^, N = 14. For N=ft=l4 the computer time was about 
four minuets and N=26, M=48 required twenty five minuets, 
(necessary for C < o.l). 
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PROGRAM 200 will solve the following problems illus- 
trated below. Note, if C = 0, A or ASTAR may not equal zero. 




•• 


— — 

- — — 

A 



C = 0, A 7^0, B?^0, ASTAR ^0, 
BSTAR^O, SIGMA0 = 0 


A?^0, Bt^O, C?^0, ASTAR 7^0, 
BSTAR 7^ 0 


A=0, Bt^O, C?^0, ASTAR 7^0, 
BSTAR 7^ 0 


A 7^0, Bt^O, Ct^O, ASTAR =0, 
BSTAR 7^ 0 



A=0, Bt^O, Ct^O, ASTAR =0, 
BSTAR 7^ 0 






